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Abstract 

We present results for matched distributions of a range of dijet event shapes at 
hadron colliders, combining next-to-leading logarithmic (NLL) accuracy in the re- 
summation exponent, next-to-next-to leading logarithmic (NNLL) accuracy in its 
expansion and next-to-leading order (NLO) accuracy in a pure Ug expansion. This 
is the first time that such a matching has been carried out for hadronic final-state 
observables at hadron colliders. We compare our results to Monte Carlo predictions, 
with and without matching to multi-parton tree-level fixed-order calculations. These 
studies suggest that hadron-collider event shapes have significant scope for constrain- 
ing both perturbative and non-perturbative aspects of hadron-collider QCD. The 
differences between various calculational methods also highlight the limits of relying 
on simultaneous variations of renormalisation and factorisation scale in making reli- 
able estimates of uncertainties in QCD predictions. We also discuss the sensitivity 
of event shapes to the topology of multi-jet events, which are expected to appear in 
many New Physics scenarios. 
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1 Introduction 



Event shapes measure the geometrical properties of the energy flow in QCD events and, 
notably, its deviation from that expected based on pure lowest order partonic predictions. 
Event shapes, as well as being among the first observables proposed to test QCD |T] , have 
been inextricably tied with the progress of QCD. They have played a crucial role in the 
extraction of the strong coupling from properties of the final-state [2]. They have been 
essential in tuning the parton showers and non-perturbative components of Monte Carlo 
event generators [3, IH El E] and have also provided a laboratory for developing and testing 
analytical insight into the hadronisation process (e.g. refs. [7] and the reviews [8l|9]). From 
a technical point of view, the development of resummations and fixed-order calculations 
has benefited from comparisons of predictions for event-shape distributions obtained with 
both kinds of methods [101 El 1121 HSl EH HSl HB]. Additionally, they are one of the several 
tools that are used for classifying hadronic final states in new physics searches. 

The majority of investigations of event shapes has been performed for e~^e~ colliders, 
with significant work also in DIS. A review of some of that work is given in [17]. In contrast, 
few dedicated studies have been performed on them at hadron colliders, with a handful of 
measurements at the Tevatron [TH [191 [20] , ^ pure fixed order study [21], pure resummations 
in [221 [231 [21], and a recent experimental simulation study by CMS [2S], as well as some 
investigation of the use of event shapes applied to jet contents for the identification of 
hadronic decays of boosted massive particles [23 [23 [21] (other approaches are reviewed 
e.g. in ref. [29]). The purpose of this article is to help bring our understanding of hadron- 
coUider event-shape phenomenology closer to the level of sophistication that is standard 
in the e~^e~ and DIS cases, concentrating specifically on event shapes in hard QCD (dijet) 
events. 

As is well known from the e^e~ and DIS cases, accurate studies of event shapes in- 
volve the simultaneous use of two kinds of calculation. Fixed-order calculations provide 
expansions of event-shape distributions in powers of the strong coupling, «<,. They are 
available up to next-leading-order (NLO) for hadron-collider event shapes, through the 
NLOJET++ [21] program. When the event shape has a value f <^ 1, for each power of 
in the distribution there can be up two powers of a large logarithm, L = Inl/v, associated 
with soft and collinear enhancements. This compromises the convergence of the pertur- 
bative series. The enhanced terms can however be resummed to all orders, providing the 
dominant contribution to the distribution for f <^ 1. Such resummed predictions tend to be 
carried out for the distribution integrated up to some value v, which generally has an expo- 
nentiated structure exp{Lgi{asL) + g2{asL) + aggsiagL) + . . The Lgi{asL) term gives 
leading logarithmic (LL) accuracy in the exponent, g2{c(sL) is next-to- leading- logarithmic 
(NLL), etc. For suitable (recursively infrared and collinear safe, global) observables, the 
CAESAR program ^31j calculates both Lgi{asL) and (72 («s-^)- To obtain a reliable predic- 
tion for an event shape distribution, one must combine both types of calculations, via a 
"matching" procedure [32]- An appropriately performed matching of NLO fixed order and 

^Not all event-shape distributions exponentiate, see e.g. jSOllSl]- 
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NLL exponentiated resummation allows one to ensure that in the expansion of the resum- 
mation one correctly accounts for all terms O-^U'-, with 2n — 2 < p < 2n, which is NNLL in 
the expansion. While NLO+NLL with NNLL in the expansion is the state-of-the-art for 
generic e"'"e~ and DIS event shapes|^ matching of this kind had not so far been achieved 
for hadron-collider event shapes. 

Our study here will give NLL-I-NLO predictions for event shape distributions both at 
Tevatron {pp, y/s = 1.96 TeV) and at the LHC {pp, y/s = 14 TeV). We start in Sec. E] 
by recalling the definitions of three classes of global event shapes for hadron collidersjfl 
addressing also the question of event shapes defined in terms of jets rather than particles. 
In Sec. [3] we describe the structure of the perturbative resummation as well as its match- 
ing to the NLO result. We also discuss possible general event-shape resummation issues 
associated with "super-leading" logarithms [391 HQ]. Sec. IHwe present our results for 
matched distributions, paying particular attention to the issue of theoretical uncertainties. 
We also compare our results to those obtained with parton-level shower Monte Carlo event 
generators, in some cases also matched to exact multi-parton tree-level matrix elements. 
In Sec. [5] we briefly discuss the impact of non-perturbative effects, the hadronisation and 
the underlying event. Finally, switching to more phenomenological questions, in Sec. |6] 
we compare various event shapes' ability to distinguish characteristically different event 
topologies, and examine their robustness in such tasks, both with respect to parton show- 
ering and to the orientation of the final state event. Our results are summarised in Sec. [71 
Some technical details are collected in Appendices [Aland El Many further additional plots 
can be obtained from the URL [41] . 



2 Event-shape definitions at hadron coUiders 

In this article, we shall consider observables that measure the extent to which an event's 
energy flow departs from a dijet structure. The lowest-order contribution to a dijet event 
consists of just two incoming and two outgoing partons. Throughout the paper we refer to 
these QCD configurations as the "Born limit" . Many of the event shapes studied here were 
presented for the first time in [23]- AH share the property of continuous globalness [HI 
[31] . which is a necessary condition for being able to carry out a resummation to NLL 
accuracy without a leading- A^^c approximation, and which also contributes to the simplicity 
of Caesar's generalised resummation approach (independently of the question of large- 

^ For specific observables, higher logarithmic and/or fixed-order accuracies have been reached. This is 
e.g. the case with NNLL accuracy in the exponent for the e+e" energy-energy correlation |3^, NNNLL 
for the thrust distribution in e+e~ |15) and NNLL for the Higgs or vector boson transverse momentum 
spectrum at hadron colliders [34l|35]. Additionally NNLO accuracy has been achieved for a range of e^e~ 
event shapes [H], with NNLO-^NLL matching in [35] and NNLO-^NNNLL in [T5] . 

For non-global observables [TT], those sensitive to emissions only in a restricted phase space region, 
the angular-ordered branching underlying CAESAr's resummations does not account for all NLL effects. 
Additional soft, large-angle contributions have to be resummed, which is currently possible only in the 
large- iVc limit [11] [37] (though see also [38j). This is the reason why, in this work, we consider only the 
global case. 
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Nc approximations). For an observable to be continuously global, it has to be sensitive 
to all emissions in an event (this is the requirement of globalness), and moreover it should 
have definite scaling properties with respect to secondary emission's transverse momenta 
(see sec. I3.1.1l for a mathematical formulation). The continuously global event shapes we 
propose fall into three main classes: observables that are directly global, others that are 
supplemented with "exponentially suppressed forward terms" and observables with "recoil 
terms" . 



2.1 Directly global observables 

We first consider observables that are defined in terms of all hadrons in the event, therefore 
the name 'directly' global. The global transverse thrust is defined as 

rr — I^J-* ■ -I ^ 

T^^g = niax — , (2.1^ 



where the sum runs over all particles qi in the final state, q±i represents the two momentum 
components transverse to the beam, q±i its modulus, and ht is the transverse vector 
that maximises the sum. The observable which is resummed is then T±^g = 1 — T±^g, 
which vanishes in the Born limit. The normalization of event shape observables to a 
hard transverse scale of the event is important because it reduces uncertainties associated 
with the experimental jet-energy scale, which partially cancel between numerator and 
denominator [l3]. For most event shapes (except r^ ^) the choice of specific hard scale to 
which one normalises is arbitrary, and could for example also be the sum of the transverse 
momenta of the two hardest jets. 

The transverse thrust axis and the beam form the so-called event plane. One can 
then define a directly global thrust minor, which is a measure of the out-of-event-plane 
energy fiow 

T^^,^^^^^. (2.2) 
In close analogy with the e~^e~ case [B], one can formulate a transverse spherocity: 

^±,, = -T . ."^in (2.3) 



■'^ 4 n=(n^,nyfl) \ ^^q^i y 

where the minimisation is carried over all possible unit transverse 2-vectors n This 
variable ranges from for pencil-like events, to a maximum of 1 for circularly symmetric 
events. 

An alternative observable, which makes use of a linearised version of the transverse 
momentum tensor (with direct analogy to the C and D parameters |35] used in e^e~), is 



"^Numerically, the minimisation is simplified by the observation (based on extensive numerical tests) 
that the n that provides the minimal sum always coincides with the transverse direction of one of the qi. 
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the i^'-parameter: 




where Ai > A2 are the two eigenvalues of M . Related variables have been considered in 
the plane transverse to the thrust axis in e^e~ (resummed for 3-jet events in |16]) and in 
the plane transverse to a jet in the context of boosted top-quark decays [261 EH ET], where 
forms involving the determinant of M^*", e.g. 4AiA2/(Ai + A2)^ = 4F/(1 + F)^, have been 
used. There is a one-to-one mapping between different forms, and we have chosen eq. (12. 4 p 
because it gives clearer separation between different kinematic regions. 

Finally, we consider the exclusive variant of the fcj-algorithm [17] (closely related to the 
inclusive variant [48] as adopted for Run II of the Tevatron [49] and expected to be used 
also at the LHC) 

1. One defines, for all n final-state (pseudo)particles still in the event, 

dkB = qlk , (2.5) 
and for each pair of final state particles 

rl ^2 1 iVk - yi? + (0fc - 00' ...X 

where ?/j = | In ^^p", is the rapidity of particle i and 0j its azimuthal angle. The 
jet-radius parameter R sets the angular reach of the jet algorithm. Throughout this 
paper, we will take i? = 0.7. 

2. One determines the minimum over k and / of the dki and the dks and calls it 

If the smallest value is diB then particle qi is included in the beam and eliminated 
from the final state particles. If the smallest value is dij then particles qi and qj 
are recombined into a pseudoparticle (jet). A number of recombination procedures 
exist. We adopt the E-scheme, in which the particle four- momenta are simply added 
together, 

qij =qi + qj ■ (2.7) 

3. The procedure is repeated until only 3 pseudoparticles are left in the final state. 
The observable we resum is the directly global three-jet resolution parameter 

?/23 = ^max{ci(")}, (2.8) 

-i£ n>3 

where P± is defined by further clustering the event until only two jets remain and 
taking Pi as the sum of the two jet transverse momenta, 

P±=P±,i+P±,2. (2.9) 
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While directly global event-shapes are defined in terms of all particles in the event, 
experimental measurements can be carried out only up to some given pseudorapidity r^max 
iVmax ~ 3.5 at the Tevatron and r/max ~ 5 at the LHC). However, as long as the event- 
shape's value V is not too small [22j, v > fmin, one can safely neglect the contribution 
of hadrons beyond the rapidity cut. The value fmin up to which a NLL resummation is 
valid is observable specific. In particular, it depends on the behaviour of each event shape 
under a soft and collinear emission, as derived in |22l |23]. Further discussion is given in 
appendix [Bl 



2.2 Observables with exponentially suppressed forward terms 

One way to address the difficulty in performing measurements near the beam is to define 
event-shapes using only particles in a central region and to add a term sensitive to emissions 
along the beam direction, so as to render them global, but with an exponential suppression 
in the forward or backward directions. We define the central region C by requiring that the 
rapidity of particles in C satisfies |?7j| < rjc = ?/j,max + ^^7, where |/j,max specifies the rapidity 
region in which the two highest pt jets should lie, and 61] is a rapidity buffer around the 
jets of size ~ 1. 

Given the central region C, we introduce the mean transverse-energy weighted rapidity 
Tjc of this region, 

^c = -Q — ^Vi(l±i, Q±,c = '^q±i, (2.10) 

^-^''^ iec i&c 

and define the exponentially suppressed (boost-invariant) forward term as 

^c—T^E^^^^"'"'"'''- (2.11) 

We can then define non-global variants of the event-shapes defined in sec. I2.1l by restricting 
the sums to just the central region. For example, we have a central transverse thrust, 

T^,c = max — ^— , rj_,c = 1 - 7^i,c , (2.12) 

a central thrust minor 

Tm,C 

iec 

and a central three-jet resolution threshold, 1/23^ defined by the algorithm of sec- 
tion [2]T] applied only to the final state particles in C (but maintaining the "beam" distance, 
eq. (123|)). 

Finally, we define "exponentially suppressed" variants of the event-shapes of sec. 12.11 
by adding to the central event-shapes a power of Sq which makes the event-shape continu- 
ously global [23]. We obtain the exponentially suppressed transverse thrust, thrust 
minor and three-jet resolution, 

r±,£ = T^,c + Sc, (2.14) 



(2.13) 

V±,c — 
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Tm,£ — Tm,C + 5 (2-15) 

y2z,£ = y2^fi + ■ (2.16) 

Additionally, one can consider event-shapes which are more naturally defined using 
particles only in a restricted region, like jet-masses and broadenings. Given a central 
transverse thrust axis nr,c, one can separate the central region C into an up part Cu 
consisting of all particles in C with pj_ ■ fiTfi > and a down part Cd with p± ■ ut.c < 
respectively. One then defines, in analogy with e~^e~ [50], the normalised squared invariant 
masses of the two regions 

V X = U,D, (2.17) 

yi€Cx / 

from which one can obtain a (non-global) central sum of masses and heavy-mass, 

Ps,c = Pu,c + PD,C , pH,c = ^ax{pu,c, pD,c} , (2.18) 

and the corresponding global event-shapes, the exponentially-suppressed sum of masses 
and heavy-mass 

ps,s = ps,c + £c 5 pH,e = pH,c + £c ■ (2.19) 

With the same division into up and down regions as for the jet masses, one can define 
jet broadenings. To do so in a boost-invariant manner, one first introduces rapidities and 
azimuthal angles of axes for the up and down regions, 

r/x,c-%^^^, ^x,c^^f^, X = U,D, (2.20) 

and defines broadenings for the two regions, 

Bx,c = Yl ^^V^''^ ~ + ('^^ ~ ' X = U,D, (2.21) 

from which one can obtain central total and wide-jet broadenings, 

Bt,c = Bjj^c + -Bd,c 5 Bwfi = max{Bu^c^ Bd c} ■ (2.22) 

Adding the forward term one obtains the global exponentially-suppressed total and 
wide-jet broadenings, 

Bt,£ = -Bt,c + £c 5 Bw,£ = Bw,c + £c ■ (2.23) 

We note that an observable that effectively has exponentially suppressed forward behaviour 
has also been studied in [241. 
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2.3 Observables with recoil term (indirectly global observables) 



Because of transverse momentum conservation, if radiation is emitted in the forward region 
C, recoil effects will cause the vector sum of the transverse momenta in the complementary, 
central region C to be non-vanishing. It is then possible to exploit this effect to make 
observables (continuously) global, despite the fact that only a central subset of particles 
in the event effectively enters the definition of the event shapes. To do so, we add to the 
central event-shapes a suitable power of a recoil term, the two-dimensional vector sum of 
the transverse momenta in C, 

(2.24) 

We obtain than recoil enhanced transverse thrust, thrust minor, three-jet 
resolution, sum- and heavy-jet masses, total and wide broadenings 



r±,n = r^,c + 7^±,c , (2.25a) 

Tm,n = T„,,c + 7^x,c , (2.25b) 

y23,n = Z/23,c + ^Ic ' (2.25c) 

P5,7^ = ps,c + T^±,c , pH,^ = pH,c + T^Lfi , (2.25d) 

Bt^-r. = Brfi + 7^±,c , -Bw,7e = By/fi + "^±,0 • (2.25e) 



2.4 Particles versus jets as inputs 

The event shapes discussed so far have all been defined in terms of the particles in the 
event. The experiments don't measure particles directly. They do, however, have methods 
such as the combination of information from electromagnetic and hadronic calorimeters 
into "Topoclusters" (ATLAS j5^) and, with tracking, "particle flow" (CMS |52]), that 
provide inputs to jet algorithms that are quite close to particles. These same inputs would 
probably also be well suited to event-shape studies. 

In uses of event shapes to cut on event topology in beyond-standard-model searches, 
as well as in the study of ref. [25], it is not particles but instead jets that have been used 
as inputs. The jets are usually defined through an angular resolution parameter R (as 
in eq. (12. 6p ) and a transverse momentum cutoff, which we will denote pto- One of the 
interests of using jets is that the pto cutoff eliminates much of the contamination from the 
underlying event, which can easily contribute O (100 GeV) of transverse momentum to the 
rapidity region covered by LHC detectors. 

From the point of view of resummation, the use of jets as inputs poses two main 
problems. One comes from the presence of the new scale pto in the problem: in terms 
of the parameters a and 61^2 defining the event-shape's sensitivity to radiation along the 
incoming legs (cf. table [T] and section I3.1.ip . this new scale causes separate regions of 
event shape value to each involve different logarithmic structure: for cases with 61 2 > 
the potentially different regions are v ^ {Vto/QYi ipm/QY ^ v ^ ijPm/QY^^^''^ and 



7^ 



9 



V ^ {pm/ QY^^^ '^ ■ The first of tliese regions may be witliin tlie scope of CAESAR if pta/Q 
is sufficiently small. 

A second problem is that of globalness. Emissions coUinear to any outgoing hard parton 
will be clustered together with its emitter to form a jet. The observable's sensitivity to 
these emissions will then depend on the jet recombination scheme. In the i?-scheme, the 
current default at the Tevatron and LHC, the jet four-momentum is constructed by simply 
adding the four momenta of its constituents. Therefore, all observables defined using 
transverse momenta will get no sensitivity to emissions inside each of the two hard jets, 
and will therefore be non-global. This statement is true for any recombination scheme 
that adds three-momenta vectorially. For variables with sensitivity to longitudinal degrees 
of freedom, globalness can only be assessed on a case by case basis. For instance if one 
considers any global version of the total and heavy-jet mass (with exponentially suppressed 
or recoil term), in the i5-scheme the mass of each central hard jet will enter the hemisphere 
central jet mass pxfi in gQ- f l^-17p . Therefore one obtains the same result for the central 
component of the event shape as would have been obtained using hadrons as inputs (modulo 
the fact that the jet clustering may affect which particles are considered central). 

One alternative to the use of jets as inputs, in order to avoid the globalness issue, is the 
following: use as inputs the particles that are inside the two hardest jets, together with 
all the remaining jet momenta. Note that this does not eliminate the issue of the extra 
scales related to pto, though it does maintain the reduced sensitivity to underlying event 
that comes from the use of jets. 



3 Structure of the perturbative calculation 

Typically one wishes to consider event shapes only for events that are sufficiently hard, 
requiring for example at least one jet above some minimum transverse momentum threshold 
Pt,min and in some central rapidity region. We will denote this kind of hardness selection 
cut by a function l-i^qi, . . . , qN) of the particles in the event; ^^(gi, . . . , qN) is equal to 
1 for events that pass the cuts and otherwise. One can then define the cross section for 
events that pass the cuts, 

a = Y, I d^^^n{q^,...,q^), (3.1) 

where daN/d^N is the differential cross section for producing particles in some config- 
uration One can determine a perturbatively as long as H corresponds to an infrared 
and collinear (IRC) safe selection procedure. 

One also defines the partial integrated cross section S(t') for events that pass the cut 
and for which additionally the event shape observable V{qi, . . . , q^) is smaller than some 
value V, 

e('" - y{<li, to)) ...,?«) . (3,2) 



N 



d^ 



N 
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The differential normalised distribution for the event shape is then given by 

1 di:{v) 

a dv 

Perturbatively we will write a and S(t') as expansions in the number of powers of the 
coupling that they contain, 

a = (To + o"! + 0-2 + • • • (3.4) 

where (Tq is the leading order (LO) result, cxi the NLO result, etc.; (Xj is proportional to 
ttg"*"*. We have chosen not to extract the powers of from the a, coefficients, because the 
scale of Us may depend on the kinematics of the events over which one has integrated. 
The expansion for S(f) is similar 

^v) = So(t;) + Si(t;) + ^^{v) + . . . , (3.5) 

with the property that So(f ) = ctq because the observable vanishes at O {oil), ^il"^) looks 
like a NLO term in eq. (13. 5p . but it is usually determined from the LO term for the 
differential cross section of 

Ei(t;) = (Ti + Ei(t;), f.,{v) = - j^dv' . (3.6) 

The quantity Ti2{v) is similarly determined from the NLO term of the differential cross 
section of v. In the following we shall never use explicitly E2, since (T2, the NNLO correction 
to the dijet cross section, has yet to be calculated and since its effect would lead to terms 
that are beyond our accuracy in differential distributions. 

3 . 1 Resummat ion 

Resummations are relevant in the region of small v, where logarithmically enhanced con- 
tributions of soft and collinear origin, as large as (a^ln^f)", appear at all orders in the 
integrated cross section S(f), thus making fixed-order predictions unreliable. There is a 
large class of observables for which one can write a common "master" resummation for- 
mula, as was done in |31], in order to sum such terms to all orders in a^. In this section 
we will first examine what the class of observables is, and then review the broad structure 
of the resummation. 

3.1.1 Prerequisites for resummation with caesar 

In order for an observable to be resummed within the CAESAR framework, its functional 
behaviour in the presence of an arbitrary number of soft and/or collinear emissions has 
to satisfy a number of conditions. These have been extensively discussed in [31], and are 
checked automatically by caesar given a computer subroutine that computes the value 
of the an observable given a set of four-momenta. The conditions are: 

11 



(3.3) 



1. a specific functional form for the observable's dependence V"({p}, k) on the momen- 
tum of a single soft emission k, colhnear to one of the hard "Born" partons ("legs") 
in the event: 

Vm,k) = dA-^\ e-'^'^'%e{<P) , (3.7) 

where {p} denote the Born momenta (including recoil effects) and k is the soft 
colhnear emission; and t]*-^-* denote respectively its transverse momentum and 
rapidity, as measured with respect to the Born parton ('leg') labelled i; (f) is the 
azimuthal angle of the emission with respect to a suitably defined event plane (when 
relevant); (/((/)) can be any function for which J d(j) In g{(f)) is well defined; Q is a hard 
scale of the problem (taken here to be the sum of the transverse momenta of the two 
hardest jets). 

2. continuous glohalness [HI |53], a requirement on the observable's single-emission scal- 
ing properties in every region of the phase space. First, all the have to be equal, 
ai = 02 = . . . = a, and the di have to be all non-zero. Second, the observable's scal- 
ing at the boundaries of the soft colhnear region has to be consistent with eq. (13. 7p . 
i.e. in the soft large-angle region we require V{{p},k) ~ k^ for a fixed angle of fc, 
whilst for hard emission colhnear to leg ^ we must have V{{p},k) ~ k^'^^' at fixed 
energy for k. 

3. recursive infrared and collinear safety, a subtle mathematical condition (see [31] for 
its precise formulation) concerning the observable's scaling in the presence of multiple 
soft /collinear emissions. 

Table [1] summarises the values of the coefficients ai and hi for the event shapes presented in 
section O We stress that central observables, like the central transverse thrust eq. (I2.12p . 
defined using only hadron momenta in a selected rapidity interval, tend to have di^2 = 0, 
and therefore be non-global, [but not true for T±^g.] The exponentially-suppressed term 
Sc in eq. (12. lip or the recoil term TZ±^c in eq. (12.240 . as explained in sections 12.21 and |2.3[ 
are added to central event shapes precisely so as to make them global. The different 
powers of Sc and 7l±,c that appear in the definition of these modified event shapes (see 
for instance eqs. (I2.14p and (12.161) ) are chosen so as to ensure their continuous globalness. 
This can be seen by observing that, for each event shape, the coefficients ai corresponding 
to different legs are equal. This above discussion holds for most observables but there 
may be exceptions. For example, the central variant of the thrust-minor Tm,c is actually a 
global observable because of an indirect sensitivity to non-central emissions due to recoil. 
This is the reason why the 61^2 coefficients for Tm^s are not those that usually appear for 
observables with exponentially-suppressed components, but are rather those typical of a 
(linear) recoil term. 

Recursive infrared and collinear (rlRC) safety, a detailed discussion of which is beyond 
the scope of the present paper, is trivially satisfied for all observables that we discuss here. 
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ai,2 &i,2 


03,4 &3,4 




1 
1 1 

1 1 


1 1 
1 1 

1 


T 

T -77 

-'- m,7<. 


1 
1 


1 n 

1 
1 


y23 

y23,£ 
1/23,11 


2 
2 2 
2 


2 
2 
2 


Bt/W,£ 

BT/w,n 


1 1 
1 


1 
1 


PS/H,£ 

ps/H,n 


1 1 
1 


1 1 
1 1 


Fg 


1 


1 1 


Qphero 


2 


2 



Table 1: Table of event shapes being considered here and the powers of their parametric 
sensitivity to the transverse momentum (a^) and collinear angle (bi) of an emission along 
incoming (ai,2,&i,2) and outgoing (03,4,63,4) hard partons. 



3.1.2 NLL resummation structure 

For global observables, in events with t> ^ 1, it is possible, unambiguously, to associate the 
event kinematics with that of a 2 —t- 2 (Born) event. This is because the requirement v 
forces all radiation to be either soft or collinear. At perturbative level it is also possible to 
unambiguously attribute a partonic subprocess to the event, for example qq — )■ qq (doing 
so requires a flavour infrared and collinear safe procedure, as in [53], but the result is 
independent of the choice of the procedure). Here we will use B to label the event's 2 — )■ 2 
kinematics and 6 to label its 2 —t- 2 flavour structure. Then we can write S(f) as a sum 
over partonic subprocesses and an integral over Born configurations that pass the hard 
event cuts, 

dB^^n{B), (t;«l). (3.8) 

o 

The ambiguities in such a decomposition of S(t>) are suppressed by powers of v. 

For observables that satisfy the properties of the previous section, the result of ref. j31] 
is that we can write 

f^^ = f!^/f(.„)(l + OK)) (3.9) 

where da^^^"* (v) / dB is the LO cross section, differential in the Born configuration, separated 
into subprocesses, and understood to have been evaluated with a factorisation scale /i^ ~ 
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Q. The function f^\v) encodes the resummation, and has the form [551 132] 



/i 



exp 



L = In 



(3.10) 



where = Q;s(/Xij), with iir some renormahsation scale of order Q. 



The order-by-order expansion of f^\v) involves terms of the form a^L'^^. It is because 
of the property of "exponentiation" (a consequence of rlRC safety and of coherenc^) that 
one can write it in the form eq. flS.lOp . whose leading- logarithmic (LL) contribution in 
the exponent, Lgf\asL), resums terms a^L""*"^: i.e. corrections to the first order a^L^ 
term involve only powers of asL. The function g2%{asL, hr, jj,p) resums "next-to-leading 
logarithmic" (NLL) terms in the exponent, a"L"', also referred to sometimes as single- 
logarithmic terms. 

The LL function Lgf\asL) can be computed analytically given only the ai and hi 
values. It is given by 



Lg^^\asL) 



? 27r/3oA6, 



[a - 2A) In 1 



[a^hi- 2A) In 1 - 



2A 



a^ht 



(<5) 



OL.U 



a(a + hi) vr 



(3.11) 



where Cf^ is the colour charge {Cp or C a) of hard parton (. for the hard-scattering subpro- 
cess 5, A = asf3oL and /3o = (llC*yi — 4TRnj)/ (127r). Since the coefficients a = ai = a2 = ■ ■ 
and hi do not depend on the particular momentum configuration B of the hard partons 
gf\asL) is also independent of B. Its dependence on the subprocess arises only througt 
the colour charges of the incoming and outgoing partons. 

The NLL function g2^^{asL, hji, fip) can be decomposed into three types of terms 



gi%{asL, ^iR, Hp) = g^^J,^{asL, fiR) + ^\n 



qf\xf\v^fip) 



qf\xf\np) 



+lnJ^^^\R'{a,L)) , (3.12) 



The term g2s{oisL) accounts for NLL corrections associated with the event kinematics, the 
particular values of the di and gi{(t)) coefficients in eq. (13. 7p . the choice of renormahsation 
scale and scheme used in in gi{asL), as well as the non-trivial colour evolution of 
large-angle soft virtual gluon resummation [561 EH EH Ell [60] . 

The second term on the right-hand side (RHS) of eq. (13.120 involves parton distribu- 
tion functions (PDF) for the parton flavours in the initial state of the given subprocess, 
/ii?), at a longitudinal momentum fraction for each leg that depends on the 
This term arises because the PDFs in da^^^ /dB in eq. (13. 9 p were eval- 



X c 



Born kinematics, 
uated at a factorisation scale ^p 



Q. The presence of a PDF at scale Hp ^ Q implies 



^The validity of coherence is brought into question by the findings of refs. 
possible implications of this in section 13.31 



I SD] and we discuss the 
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that one integrates over all possible incoming coUinear emissions, up to kt ~ Q. However 
the requirement that the event shape be small, V{k) < v, restricts collinear emissions to 
have kt < v^^^"'~^''^^Q. Thus the PDFs should actually be evaluated at a factorisation scale 
~ ~ v^^^"'^'''^ hf (as occurs also in Drell-Yan transverse-momentum resumma- 

tions EB E2])- The ratio of PDFs in eq. f l3.12p serves to replace qf \x^f \nF) as used 
in the Born cross section with a PDF at the correct factorisation scale. 

The third term on the RHS of eq. (13.121) accounts for the NLL corrections associated 
with the presence of multiple soft and collinear emissions, when each has V{k) ~ v and 
they are all widely separated in rapidity. It is a function of 

R\asL) = -dLLgi{a,L), (3.13) 

and is known analytically for some observables (e.g. T±^g), while in all other cases CAESAR 
can compute it numerically via a suitable Monte Carlo procedure. J-'^^\R') sometimes 
depends on the underlying scattering channel 6, but not (for the observables studied here) 
on the hard momentum configuration. 

The behaviour of J^{R') with increasing R' (decreasing v) is a characteristic feature of 
each event shape. It depends on whether multiple emissions tend to increase or decrease 
the value of the event shape. In the first case, for a fixed value v, J-'{R') has to account 
for an extra suppression of emissions so as to keep the event shape's value less than v, i.e. 
J'iR') < 1. For the special case of V{ki, /c2, • • •) = max(l^(fci), V{k2), ■ ■ ■) then J'(-R') = 1 
(for example the 7/23 jet resolution threshold [12] in the e'^e" Cambridge jet algorithm [63]). 
Conversely if the contributions of multiple emissions tend to cancel, the function J^{R') 
has to compensate the excessive suppression given by the LL function Lgi{asL), therefore 
J^{R') > 1. 

This last case appears most dramatically when it is a cancellation between multiple 
emissions, and not a direct veto on real emissions, that is the dominant effect that keeps the 
event shape small. In this case the LL function Lgi{asL) (whose functional form depends 
only on the effect of single emission) no longer accounts for the dominant contribution to 
the distribution. Furthermore, no NLL function such as J^{R') can fully compensate for 
this. This inconsistency reveals itself through a divergence of J^{R') at a given critical value 
R'^, which can be inferred from considerations on the cancellation mechanism, as explained 
in refs. [23l[31]. Such a divergence is present, for example, for Tm,c and T^^e and occurs at 
R'^ = Ct/ {Ci + C2), where Ci and C2 are the colour charges of incoming partons and Ct 
the total colour charge of the hard parton system. It will prevent us from obtaining sensible 
NLL resummed results for these observables. In the case of recoil observables there is also 
a divergence, but at larger R'^ (smaller v), for example at R'^ = 2Ct/{Ci + C2) for Tm,n 
and Bt/w^tz- The effect on the corresponding differential distributions will be discussed 
when presenting matched results. 

3.1.3 NNLLs accuracy 

As well as discussing the LL, NLL, etc. accuracy of resummation in the exponent of 
eq. fl3.10p . one can also discuss the accuracy in the order by order expansion of S itself. 
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In this way of counting logarithms, "LLs" terms involve powers NLLs involve 

q;"L^"~^, etc. A NLL resummation in the exponent automatically guarantees NLL^ ac- 
curacy. However it is also possible (if not entirely straightforward), given the information 
at our disposal, to obtain NNLL^ accuracy. To see how, observe that the terms that we 
neglect in eqs. (13. 9p and (I3.10p are an overall as correction without logarithms, as well as 
terms a^L""^ in the exponent, starting at a^L. The latter, if they multiply a^L^" when 
expanding the exponent, lead at most to a^^'^L?"'^^ ~ a^L^""^, i.e. they are NNNLL^. 
We can therefore ignore them. As for the overall correction, when multiplied by the 
double logarithms, it gives us terms ttg+^L^" ~ a^L^""^, which are NNLLs and therefore 
cannot be neglected. This means that we need to determine the coefficient of the pure 
O {as) term. 

To do so, let us define the NLL resummed cross section as 



dB dB 

'NLL,/3(d) 



"^-"^I^Lsiv). (3.14) 



with f^}^ „, containing only the LL and NLL resummation terms. 



/nll,b(^) = exp \ Lg'(\asL) + g)°l^{asL, ^r, ^f)\ ■ (3.15) 

Then we can determine the coefficient Cf\ in terms of the first order expansion of the 
exact and resummed distributions, and in particular their difference as — )■ 0, 

„.C-.„.(^-i^)/=f . (3,16) 
v^o \ dB dB j dB 

The cf\ constant involves many contributions, including parts that cancel the and 
dependence present in the Born cross section, parts that are sensitive the observable's 
exact behaviour with respect to soft large-angle emission and hard coUinear splitting and 
parts related to the exact structure of the 1-loop 2 — )■ 2 scattering diagram. 
Now we can write the NNLL^ resummed distribution as 

(1 + asC^^o) . (3.17) 

The fact that we may multiply (1 + a^C^'^g) and '^^^g''"'' in order to get NNLL^ accuracy 
is a consequence of the property that soft-coUinear virtual corrections, which give powers 
of agL'^i affect neither the flavour, the momentum nor the colour involved in the hard 
scattering or the PDFs and therefore straightforwardly multiply all the more complicated 
contributions that are present in Cf^. 

In practice, it is not feasible to (numerically) determine the exact flrst order distribution 
for V fully differentially in the Born conflgurations. Considering instead quantities S*^^)(f), 
Er^^(f) (and their order- by-order expansions), integrated over conflgurations that pass the 
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event cuts, as in eq. f l3.8p . one can define a Ci"^ coefficient averaged over Born momentum 
configurations, 



{asCl >) = hm = ^ / dB-—asC\ >, . (3.18) 

Writing 



(l + (a,Cf)))S(^)(t;), (3.19) 

gives a distribution that is still correct to NNLL^ accuracy, because the LL, e'x.]i{Lgf'' (agL)), 
component of Tii^\v) is independent of the momentum configuration. 

In contrast, if one considers Ci averaged additionally over subprocesses 

{asCi} = hm = — > iasC{ ')cr^ , (3.20) 

then 

(l + (a,Ci))S,(t;), (3.21) 

is not accurate to NNLL^, because 

X^(a,CfVS'^exp(L^P(«,L)) ^ (^^ ^("^^i'VS'^ j (^^^ aj^^ exp(L^f )(a,L)) 

(3.22) 

since the coefficient of the double logarithms in Lgf^ {ctsL) does depend on the subprocess, 
through the colour charges of the hard partons. 



3.2 Matching of NLL to NLO 

While the resummation of logarithms is necessary in the region where event-shape values 
are small and their logarithms large, the region of large values of V is dominated by events 
with three or more well separated jets. Those types of events are described more reliably by 
fixed order calculations. It has therefore become standard to match resummed calculations 
to next-to-leading order (NLO) to have a reliable prediction over a larger range of values 
of V. In this section we will present the formulae we use to perform the matching. 

In the following we will denote with f{v) = Tj{v)/a the integrated event-shape fraction, 
where E(t>) and a are defined in eqs. (13. 2p and (13. ip respectively. After a NLL-I-NLO 
matching this quantity should satisfy the following requirements 

1. it should respect the physical constraints that, when the event shape reaches its 
maximum value Wmax, we have /(fmax) = 1 exactly and ^^^\ _ =0; 

2. its expansion up to relative O (af ) should reproduce the exact NLO result for the 
corresponding differential distribution; 
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3. one should obtain NNLLg accuracy, i.e. all logarithms O {a^U^) with m > 2n — 2 
should be correctly accounted for, which implies that the matching formula should 
reduce to eq. (13.19^ in the limit of small v. Preferably this should be the case without 
having to go through the tedious procedure of manually determining cf^ separately 
for each event shape. 

There are various matching procedures that satisfy these requirements and therefore for- 
mally have the same accuracy. We consider here the so-called log-i? [32j and multiplica- 
tive [42] matching schemes, adapted to hadron-hadron collisions. In particular both need 
to be modified in accordance with the need, section [3.1.31 to have the O (a^) constant cf^ 
term multiply the resummation separately for each subprocess. Actually, what matters is 
not so much the subprocess but the colour charges of the incoming and the outgoing Born 
partons. Therefore we can consider all subprocesses qq — )■ gg, qq' — )■ qq' , qq — )■ qq, etc., 
with the same incoming and the same outgoing colour charges as belonging to a single 
colour channel a = gg — )■ gg. The other colour channels are qg qg, qq gg, gg — )■ gg 
and gg ^ ggE 

We denote by ^^.'^^(K) the expansion of the resummed cross section corresponding to a 

specific colour channel and by ^■"''(V") the corresponding exact fixed order prediction. In 
analogy with eq. (13. 8p . where the index a is omitted a sum over all possible colour channels 
is understood. The index i denotes the order in of the expansion (relative to the Born 
cross section). 

We have obtained fixed order cross-sections using the code NLOJET++ [21]. The publicly 
available version computes cross-sections summed over the flavour of outgoing partons. We 
therefore extended it so as to have access to the flavour of both incoming and outgoing 
partons in the calculation of ctq, (Ji and T,i{v), though not for T,2{v) since its colour- 
channel separation is not needed for NNLLs accuracy. To assign events with more than 
two outgoing partons to a definite 2 — )■ 2 colour channel we used the exclusive flavour- fc^ 
algorithm of [M] to cluster events to a 2 —i- 2 topology while keeping track of flavour in an 
infrared safe manner. During the clustering procedure, quarks of different flavour might 
end up in the same jet, giving rise to multi-flavoured jets, i.e. jets whose flavour does not 
correspond to any QCD parton. These events, which do not correspond to any Born 2 — )■ 2 
processes and have vanishing weights for v 0, will be labelled as having a = other. 

The matching equations are defined in terms of the following resummed distribution 
(and its fixed order expansions ^[."^(f ), S^"2(f)), 



^ We study here only observables whose double logarithms depend only on the total colour charge of 
the two incoming and two outgoing partons, so that we do not to distinguish incoming partons 1 and 2 
(or outgoing partons 3 and 4). This means that for the matching only the colour structure is relevant, 
therefore this colour labelling does not distinguish quarks from anti-quarks or quarks of different flavour. 
A given colour channel a is then in general a sum over multiple partonic channels 6. 




(3.23) 
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where f^\v) is the resummed probabihty f^'{v) (eq. (13.151) ) with L replaced by [321 Ii2] 



L = -ln 

P 



XyVn 



+ 1 



Xy — X ■ X\ 



(3.24) 



where fmax is the maximum kinematically allowed value of the event shape, so that L{v = 
Vraax) = 0. We take the values of Umax from the NLO calculation, which is sensible since 
we want the differential distributions to reproduce the NLO result at high v. The factors 
Xy and p modify the definition of the logarithm that one is resumming. The main effect of 
Xy is to modify the logarithm at small values of V, and it will therefore affect subleading 
logarithmic terms (the change at NLL is cancelled via a suitable compensatory term in 
g2{asL)). The main effect of p on the contrary is to modify L at large values of V , it will 
therefore mainly affect power suppressed terms. Our default values for xy, Xy are given 
by (see Appendix A of ref. ^23j) are fixed by setting X = 1 and 



(3.25) 



In the same way as one varies renormalization and factorization scales around a central 
value by a factor of 2, we will probe the xy dependence through a variation of X in the 
range 1/2 < X < 2. This will provide an estimate of the error associated with unknown 
NNLL contributions to the resummation. In principle, one can also vary the power p 
around the value 1 (as a probe of terms that are suppressed by powers of f), though for 
simplicity in the following we just fix p = 1 and therefore do not include any uncertainty 
related to its variation. 

We now introduce the log-R matching formula 



with 



y) exp 



/'(^max) 



(«)/ 



(3.26) 



(a) 



X exp 



S2(f ) - T.r,2{v) 1 



aT^other 



2a, 



(a) 



where Ti2{v) has been introduced after eq. (13. 6p . It is straightforward to verify that with 
this matching equation f{v) satisfies all three requirements listed at the beginning of this 
section. 
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The alternative, multiplicative matching (mod-R) scheme that we use is 



^ OT^other 



(a) 



1 ^ Sr(.)-^Ei?(.) 



a' pother 



CT, 



(a') 



(3.28) 



where Z = ^1 — j • This matching equation has the same matching accuracy as 

eq. f l3.26p . so that using both matching procedures provides an additional way of esti- 
mating the uncertainty in the matched distributions. 

In both matching formulae, Sr,i(f) and Sr,2(^^) require the calculation of the order as 
and expansions of the ratios of PDFs at different scales that appear in eq. (13.121) . These 
have been obtained using hoppet 

In the following we will present results for normalised differential distributions 



1 da{v) df{v) 



a dv 



dv 



(3.29) 



Notice that two- loop corrections to o"2, currently unknown, are not needed for a second 
order matching, as they do not contribute to the differential distribution within the target 
accuracy. 



3.3 Coherence- violating (super-leading) logarithms 

One of the assumptions that enters into the derivation of the generalised resummations of 
[55] is "coherence" [SB], the property that real emissions and virtual corrections at large 
angles are independent of the structure of real emissions that have occurred at small angles 
(with respect to any of the incoming and outgoing legs). Physically this can be understood 
as arising because a large-angle emission (or virtual correction) sees only the sum of colour 
charges of a bunch of coUinear partons and that sum of colour charges is conserved under 
collinear splitting^ 

The assumption of coherence is challenged by the results of ref. [5^ HO] , which found 
"super- leading logarithms" (SLL), terms that go as a^L^, when calculating the probability 
of there being no soft radiation (above scale Qe~^) in a finite patch of rapidity and azimuth. 
Based on coherence, one would have expected only terms a^U^ with m < n for such 
an observable. Therefore, one might also call the terms of [5U| WJ\ "coherence-violating 
logarithms" (CVL), a name that is suitable also in the case of observables whose leading- 
logarithmic structure involves double logarithmic terms a^L^" (for which a^L^ is not super- 
leading) . 

^For initial-state splittings, large-angle emission sees the difference in colour charges between incoming 
and outgoing partons that are collinear to an incoming direction. 
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The interpretation of the result in ref. [391 SO] is that one specific class of (soft) single 
logarithmic virtual correction, "Coulomb-gluon exchange," can be affected by small-angle 
(coUinear) initial-state gluon emission, independently of how small that angle is. This is 
because in the calculation of [391 SO] Coulomb gluons are exchanged either between two 
incoming partons or between two outgoing partons but not between one incoming and one 
outgoing parton (whereas other classes of soft contribution treat incoming and outgoing 
partons on an equal footing); real initial-state splittings, however small in angle, lead to a 
redistribution in colour between incoming and outgoing states and therefore Coulomb-gluon 
exchange cares about them (but not about the corresponding coUinear virtual initial-state 
corrections). This means that the coefficient of the Coulomb single logarithms a^L" is 
proportional to the probability of soft-collinear initial state emission, a^L^*" and hence 
one obtains terms Q,^+™-/,"-+2m^ which are super-leading with respect to the expected a^L". 

The calculation of the impact of this effect requires that one follow through the soft 
colour evolution of the 2 — )■ 2 scattering [^Sl EZl EH ESI ED], for which the Coulomb- 
exchange terms provide imaginary contributions. For the purposes of our discussion here 
it is not necessary to enter into the full detail of the soft colour evolution. Rather it suffices 
to be aware, following [ID], that for the case of vetoing emissions into a finite patch (gap) 
the lowest order coherence- violating terms are contributions with structures such as 

■ &{Q - kn)Q{kn - kt2)Q{h2 - kt3)e{kt3 - k^Mku - Qe"^) (3.30a) 
= -^atL'^ + O {atL') , (3.30b) 

where we have shown only one of the two orderings given in [lO] (the other gives either 
the same number or fewer logarithms, depending on the observable). In the integration 
measures, we have labelled each momentum with (f) if it can only be virtual, and {r/v) 
if we are considering the difference between real and virtual cases. Gluon 2, the coUinear, 
possibly real gluon, can have an angle corresponding to anywhere outside the gap region, 
down to the smallest kinematically allowed angles 62 ~ ki^lQ- Gluons 1, 3 and 4 have 
only transverse momentum integrations because they are either Coulomb exchange gluons 
or the virtual counterparts of large-angle soft-gluon emission. (In the other ordering it is 
gluon 1 that is coUinear and possibly real). The integral for gluon 4 is limited to be above 
Qe~^ because below that scale the observable places no constraint on real emissions and 
so all real and virtual effects should cancel, by virtue of unitarity. Finally, the constant C 
depends on the kinematics of the hard scattering and the definition of the gap region. 

The extension to the event-shapes case involves restricting the {r/v) integration for 
gluon 2 to regions of phase-space that are consistent with the the real gluon's contribution 
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to the event shape being < e ^. Eq. fl3.30p therefore becomes 

CVL^hp -Cat ^ ^^e(l-e2)e(Qe2-A:,2)e(Qe-^/'^-fc,2^5/'^) ^ ^ 
J kn kt2 O2 hs hi 

■ e(g - ka)eikn - kt2)eih2 - hs)eika - h,)eih, - ge"^/") , (3.31) 

where the coefficients a and b are those that appear in eq. (13. 7p for the incoming legs, 
£ = 1,2, for simphcity we have neglected the di and gi factors there, and the constant 
C may differ somewhat from that in eq. fl3.30p (since there it could depend on the gap 
definition). Three cases arise: 



b<0 


— > 


CVLcv-shp 


~ Ca^L 


(3.32a) 


b = 


— > 


CVLcv-shp 




(3.32b) 


b>0 




CVLcv-shp 




(3.32c) 



where for the cases b < the number of powers of L is that obtained when relaxing the 
constraint Qe~^l°- > /c(2^2''" to become Qe^^^'^ > kt20'^J"' (consistent with the fact that we 
have ignored factors of (i^, g^). In words, the CVL contributions only appear with a large 
number of logarithms when the collinear gluon, 2, if real, is allowed to be harder than the 
virtual Coulomb gluons (3, 4)j^ For observables with 6 < the one logarithm arises from 
the integration over kn-, while the Qe~^/°' > ^42^2''° constraint forces ^2 to be at large angles 
with kt2 ~ Qe'^l"" ^ with the knock-on effect that kfi and /ct4 should also be ~ Qe'^l"" . For 
6 = 0, we instead have just Qe^^l"" > kt2 for gluon 2, and an extra logarithm then arises 
from the integration of Q2 in the collinear region. 

The expectation for yet higher order terms is that for 6 < (all observables of this paper 
except the exponentially-suppressed ones) the results in eq. f l3.32p could be multiplied by 
additional powers of a^L giving at worst a series a"L"~^, which is subleading both with 
respect to our NLL accuracy in the exponent and to our NNLLs accuracy in its expansion]^ 
For 6 > (the exponentially suppressed observables) one would obtain terms agL^""^. It 
is not clear how they would fit into the exponential resummation, except that they would 
certainly destroy NLL accuracy in the exponent; in the expansion of the resummation they 
would represent terms NNNLL^ and therefore be subleading with respect to our accuracy. 

One caveat with regard to the above discussion is that the results of [321 SD] have been 
obtained in a strongly ordered eikonal approximation, with the assumption that the "strong 



^ Note that for the exponentiahy-suppressed observables, with 6 = a, we expect the coefficient of 
the CVL to be significantly suppressed as compared to the gap case, because of the way in which the 
event-shape constraint restricts the phase-space integration region. 

^This result involves the assumption that there must be at least two large-angle or Coulomb virtual 
gluons softer than the collinear gluon. While this is the case for the contributions found in the gap case 
[39l l40] we do not show here that it will always necessarily be the case for event-shapes. If there could 
instead be coherence-violating contributions with just one large-angle or Coulomb virtual gluon that is 
softer than the collinear one, one might expect terms up to q;"L"~^, which, however, are still subleading 
relative to the calculations of this paper. 
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ordering" parameter is transverse momentum. With different assumptions, the results 
change. For example, if the correct ordering parameter were energy, then eq. (13.301) would 
be modified in such a way as to give an infinite result. If instead one considered emission- 
time (or virtuality) ordering^ which leads to ka ^ kt202 ^ fcts, then the contribution of 
eq. (I3.30p would be halved, 

(.VL(^ordercd) ^ |^4^5 ^ (3.33) 

and the corresponding result for the event shapes case would become 

b<a ^ CVL^X^-'^ - CatL (3.34a) 
b = a ^ CVClp^"^^") ~ C-a^L^ (3.34b) 
b>a CVL^;^-^'""'^ CatL' (3.34c) 

In this case, for all the observables being discussed in this paper, the CVL terms would be 
subleading relative to our accuracy. 

To conclude: given today's knowledge it is not clear whether coherence- violating terms 
matter at our accuracy for event-shape resummations. The critical issue is that of the 
appropriate ordering parameter [kt, time or virtuality ordering, or some other ordering). 
The correct ordering needs to be derived (by going beyond the eikonal approximation), 
unless of course there exists some yet-to-be found contribution that cancels the CVL terms. 
If CVL terms do exist and kt ordering is correct, then they will invalidate our statement of 
NLL accuracy in the exponent for the exponentially suppressed class of observables, though 
not our statement of NNLL^ accuracy in the expansion of the distribution. In practice we 
have reason to believe that their numerical impact will still be small: partly because the 
CVL terms were already not very large in [32]; and partly because the large colour factors 
multiplying the double logarithms of our resummation force the majority of events to be 
in a region where the logarithms are not actually all that large (a reflection of this will 
appear in section 14.51 where we will see that naive exponentiation of the NLO calculation 
is not too different from the full NLO+NLL result, even though it misses classes of LL 
terms in the exponent and LL^ terms in its expansion). 



4 Perturbative results 

In this section, we shall consider numerical results both for Tevatron and LHC collision 
scenarios. We will start by presenting the event selection cuts that we use. We shall 
show results for NLL+NLO matched calculations for a range of observables. We will pay 

^''One gets the same result for the ordering based on two different considerations. Physically, the time 
scale for the collinear emission to take place is {l/kt2) ■ (1^2/^*2) ^ 1/(^*2^2), where lu2 is the energy of 
gluon 2, to be compared with 1/kt for a large-angle virtual gluon exchange. In terms of the diagrammatic 
structure, ordering is in part related to the virtualities of propagators and for a hard scattering of partons 
with energy E, the squared propagator virtuality induced by soft and collinear gluon emission is ^ EUJ2O2 — 
Ekt202 to be compared with Ekt for large- angle virtual gluon exchange. 
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particular attention to the estimation of uncertainties on our predictions, and comparisons 
to separate pure NLO and NLL calculations. We will also compare our results to Monte 
Carlo parton-shower results, with and without tree-level matrix element matching. 

4.1 Event selection cuts 

The Tevatron scenarios involve pp collisions at centre-of-mass energy ^/s = 1.96 TeV. 
Events are clustered with the SISCone jet algorithm [67] (similar to the MidPoint algo- 
rithm |49j that is in widespread use at the Tevatron, but infrared safe), with a jet radius 
R = 0.7 and a split-merge overlap threshold / = 0.75. The two hardest (highest-pt) jets 
in the event should have rapidities \y\ < 0.7. Events are accepted for a low-pt sample if the 
hardest jet has pt > 50 GeV, while they are accepted for a high-pj sample if the hardest jet 
has Pt > 200 GeV. As concerns the event shapes, the central region is defined by rjc = 1. 

The LHC scenarios involve pp collisions at a centre-of-mass energy y/s = U TeV0 
Events are clustered with the kt jet algorithm [HI HE], with a jet radius R = 0.7. The two 
hardest jets in the event should have rapidities \y\ < 1. Events are accepted for a low-p^ 
sample if the hardest jet has pt > 200 GeV, while they are accepted for a high-p^ sample if 
the hardest jet has pt > 1 TeV. As concerns the event shapes, the central region is defined 
by Vc = 1-5. The larger choice than at the Tevatron reflects the LHC detectors' larger 
overall rapidity coverage. 

The cross sections for the different selections are given in table [21 These, and all 
other NLO calculations presented here, have been carried out with NLOJET++ 3.0 [21] 
(modified to provide access to parton flavour information up to 0(q;^)), with CTEQ6M 
Parton Distribution Functions (PDFs) [68j and FastJet 2.3 [69] for the jet clustering. 
The renormalisation and factorisation scales have central values HR,fiF = Pt = {pti + 
Pt2)/2, where pa and pt2 are the transverse momenta of the hardest and second hardest 
jet respectively. The quoted errors correspond to the uncertainties due to scale variation 
Pt/2 < fin, fip < 2pt, with (1^/2 < fiR < 2fiF. 

Given the cross sections in table [2]one easily concludes that with available (anticipated) 
luminosities at the Tevatron (LHC) there will be large event samples on which to study 
event shapes. One also observes that NLO corrections are larger than one is used to seeing 
for (say) inclusive jet cross sections. This is a consequence of our selection based on the 
value of pti. A selection based on the average pt of the two hardest jets would instead have 
given K-factors rather similar to those for the inclusive cross section!^ 

^-'^The LHC will initially run at centre-of-mass energies that are below ^/s — 14 TeV, though the exact 
energy of collisions is subject to uncertainty and will vary over the course of the initial runs. Given that 
the generation of the NLO results for a single combination of collider energy and event-selection cuts 
requires many CPU- years of computing time, we have decided to remain with ^/s = 14 TeV as our default 
choice for the time being. The general picture as it applies to other centre-of-mass energies can be largely 
understood by interpolation between the Tevatron and 14 TeV LHC results. 

^^The choice of a cut on pti was originally motivated by the observation in the context of HERA [701 
1711 172] that identical simultaneous cuts on pti and pt2 led to poor convergence of the perturbative series, 
for reasons discussed in [731 IZll- Cutting on pti was intended as a way of avoiding this problem, but, as 
we see here, seems to introduce issues of its own. Note that it is probably not advisable to introduce a 
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Tevatron, pa > 200 GeV 
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101+i pb 


41% 


43% 


12% 


14 TeV LHC, pa > 200 GeV 


13.3+1^ nb 


23.8^11 nb 


7% 


40% 


50% 


14 TeV LHC, pa > 1 TeV 


6.4^2-0 pb 


10.511^ pb 


31% 


51% 


17% 



Table 2: Cross section for events that pass the selections cuts described in the text. The 
uncertainty is that due to scale variation with the choice pt/2 < fj,R,^F < 2pt, with 
IJ'f/2 < fiR < 2 fip, where pt is the average of the transverse momenta of the two hardest 
jets. Also shown is the breakdown (at LO) into the main scattering channels; q denotes 
both quarks and antiquarks, and channels that contribute neghgibly, such as gg qq, are 
not shown. 

Table |2] also shows the breakdown into the 3 main partonic scattering channels (as 
calculated at LO). At each of the colliders, for the lower pt cut, channels involving gluons 
are dominant, while for the higher pt cut channels involving quarks play a bigger role. This 
difference between low and high-p^ samples will be clearly visible in the final results. 

4.2 Resummed results and uncertainty studies 

Here we present resummed results for the global thrust minor, T^^g together with a study 
of its perturbative uncertainties at the Tevatron for the high-pf sample {pn > 200 GeV). 

Fig. [1^ illustrates the NLL+NLO matched distribution obtained in the log-R matching 
scheme with X = 1. The renormalisation and factorisation scales are chosen, event- 
by-event, to be fJ^F = fi'R = Pt = {pti + Pt2)/2 in both the resummation and the NLO 
calculation. The distribution has a peak at small event-shape values that is characteristic 
of all resummed (and physical) event-shape distributions. 

The uncertainty on the prediction is almost as important as the result itself, especially 
as it will allow us to gauge the significance of any disagreements that we will see with other 
predictive methods. 

The most widely used form of uncertainty estimate is the variation of renormalization 
and factorization scales. The solid (red) curves in fig. [T]d illustrate the effect of varying these 
scales simultaneously, showing the ratio of results with fip = hr = pt/2 and fip = fJ'R = 2pt 
to the default result. Except at very small event-shape values or at very large ones, 
where the distribution vanishes, one sees that the impact of symmetric scale variation 
is only about 5%. Asymmetric scale variations are shown by the dashed (green) curves, 
corresponding to fip = {pt/2,2pt} while keeping /x/j = pt, and fiR = {pt/2,2pt} while 
keeping hf = pt- For moderate and large values of the event shape they have a significantly 
larger impact than symmetric scale variations, of the order of 10% for moderate Tm,g- This 

staggered cut on pti and pt2, e.g. pti > 50 GeV and pt2 > 40 GeV, because that introduces an extra small 
parameter in the problem, related to the difference between the two pt cuts. 
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Figure 1: a) NLL+NLO resummed matched distribution for the directly global thrust 
minor at the Tevatron {pp, ^/s = 1960 GeV) with pt > 200 GeV; b) renormalisation 
and factorisation scale uncertainties for xr,xf = 0.5,1.0,2.0, varied separately, with the 
condition 0.5 < xf/xr < 2; c) effect of varying X = 0.5,0.7, 1.0, 1.5,2.0 in eq. ([321; d) 
effect of changing the matching scheme. 



26 



highlights the importance of considering both symmetric and asymmetric variations. 

Fig. [It shows the impact of varying X in eq. fl3.24p , with the hne thickness increasing 
from X = 0.5 to X = 2. As discussed in Sec. 13.21 this variation can be used to estimate 
the effect of higher order logarithms not included in our NLL resummation. We find 
that for moderate and large values of Tm,g the effect is similar in size to the asymmetric 
renormalization/factorization scale variation. Closer to the peak of the distribution (where 
the bulk of events sits) , the impact of the X-scale variation is mildly larger. We also note 
that the variation is quite asymmetric: smaller X values distort the central distribution 
much more than larger values. 

Finally, in Fig. [T]i we estimate uncertainties that arise from the details of the matching 
procedure. In particular we show the ratio of the mod-R matched distribution to the log-R 
(see eqs. (13. 26 p . (I3.27P and (I3.28P ). It is clear that at large values of T^^g, the difference 
between the two matched distributions is large, with differences of up to 45% for T^^g ~ 0.5. 
These very large discrepancies occur however only in the tail of the distribution, where few 
events are present. Comparison to NLO at high T^^g (not shown) indicates that of the 
two matching schemes, log-R matching is the one with smaller higher order terms (i.e. its 
NLO+NLL result is closer to NLO) at large Tm,g. In the following we will therefore use 
log-R matching as our default. 

The above findings are representative of the results for the other event shapes considered 
here, both at the Tevatron and at the LHC and for the low- and high-pt samples (further 
NLO-I-NLL results are shown in Sec. l4.4l and on the website associated with this article |41]). 
In particular symmetric renormalization and factorization scale variation (as is currently 
done in many phenomenological studies) systematically underestimates the true size of 
theoretical uncertainties. While detailed error-estimate studies such as variation of X- 
scale and matching procedure are possible only for specific (resummed) calculations, an 
asymmetric fiR and fip variation can be carried out for generic observables. This is perhaps 
most relevant for multi-scale observables, where the scale of as is a priori not clear. 

4.3 Comparison of resummed, NLO and matched results 

In this section we compare various levels of fixed order calculations (LO, NLO), pure 
resummed ones (NLL) and matched ones (NLO+NLL) at the Tevatron for the high-pt 
sample. Because NLL+NLO resummations are rarely available (e.g. they are currently 
not available for non-global observables), we discuss in particular the extent to which NLO 
alone can be used to describe event shape distributions. As in the previous section we will 
use Tm,g to illustrate our findings, but results are fairly independent of the specific event 
shape, the collider and the details of the hard cuts. 

Fig. |2] shows the result for the log-R matched Tm,g distribution compared to pure 
resummation (left) and pure NLO, and LO for reference (right). For the matched resummed 
result, the band corresponds to the span of all the uncertainties shown in fig. [H while in 
the fixed order calculations it corresponds to the variations of just the renormalisation 
and factorisation scales {pt/2 < Hr^Hf < 2pf, with /iF/2 < < 2 fip)', for the pure 
resummed result the band corresponds to the renormalisation, factorisation and X-scale 
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Figure 2: The distribution for the representative observable Tm,g, comparing pure resum- 
mation (NLL, left) and pure fixed order (LO and NLO, right) with the matched prediction 
(NLL+NLO). For the matched resummed result, the band corresponds to the span of all 
the uncertainties shown in fig. [H while in the fixed order calculations it corresponds to the 
(asymmetric) variations of just the renormalisation and factorisation scales; for the pure 
resummed result the band corresponds to the renormalisation, factorisation and X-scale 
uncertainties. See text for more details. 



uncertainties. 

As expected, the matched distribution agrees with the NLO results at large values of 
Tjn^g. However for the pure NLL resummation without any coefficient function obtained 
from eq. (13. 23 p . the level of agreement with NLO+NLL is quite poor even at fairly small 
values of Tm,g- For example, the position of the peak of the distribution is not all that well 
predicted (at Tm,g ~ 0.09 — 0.11 rather than at Tm,g ~ 0.08). As far as the height of the 
peak is concerned, both NLL and NLO+NLL distributions are normalized to one, however 
the NLL distribution becomes negative at Tm,g > 0.35, and this negative tail causes the 
distribution to be far too high at low Tm,g- It is on the other hand reassuring that these 
large differences with the matched distribution are reflected in the very large uncertainty 
band of the NLL distribution. 

As far as the fixed-order results are concerned. Fig. [2}d, they are as expected divergent 
at small Tm,g- The LO distribution essentially never agrees with the matched distribution, 
while the NLO does within uncertainties for Tm,g ^ 0.2. It is on the other hand evident 
that scale uncertainties of the NLO results at small Tm,g underestimate the size of higher 
order corrections not included in the fixed order calculations. 

Altogether, figure Fig.[2]highlights how neither NLO nor resummation alone can provide 
a sensible prediction, while the combination of NLO and resummation gives significantly 
reduced scale-dependence compared to either on its own. Furthermore the matching pro- 
cedure gives the general shape that is associated with the resummation, while maintaining 
the large-w behaviour of the NLO prediction. We finally note that for the event shapes 
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Figure 3: The NLO+NLL matched distribution for the directly global thrust minor, Tm,g 
shown for the Tevatron (left) and i/s = 14 TeV LHC (right) with two transverse momen- 
tum cuts for the event selection. 



presented here a LO+NLL resummation agrees well neither with the LO, nor with the 
NLO+NLL matched resuh. 

4.4 NLL+NLO matched results for a range of observables 

In the previous section we established that contrary to NLO or NLL alone, NLO+NLL 
provides robust theoretical predictions for event shapes distributions over a large range of 
the event-shape values. This section contains the bulk of results of the present work: we 
discuss NLL+NLO resummed distributions for a number of event shapes variables, both 
at the Tevatron and at the LHC and for both, low- and high-pf samples. 

We start by looking at the effect of changing the hard selection cuts and the collider 
(Tevatron vs LHC) for the same observable discussed previously, T^.g, Fig. |3l There is 
a striking similarity between the Tevatron (left) and the LHC plot (right), both for the 
low-pt (Tevatron, pti > 50 GeV and LHC, pti > 200 GeV) and for the high-p^ (Tevatron, 
Pti > 200 GeV and LHC, pu > 1000 GeV) samples. We also notice that low-pt curves 
are broader and peaked at a higher value of T^^g. This is a consequence of the higher 
prevalence of gluons in both the initial and final states of the hard scattering. 

Because of this similarly between low- and high-pt samples at the two colliders, we 
examine results for a large range of observables, as defined in Sec. [21 just for the high-pj 
cuts at the Tevatron, Fig.HJ and for the low-pt cuts at the LHC, Fig.O For each observable, 
we give two uncertainty bands, one corresponding to a symmetric scale variation (hatched, 
dark blue) and one defined in terms of all theoretical uncertainties as discussed in Sec. 14.21 
(solid, light blue). Comparing Fig. HJand Fig. [5] we see that, as observed earlier for Tm,g, 
the peaks of the distributions are further to the right and the distributions are broader 
for the LHC (low pt) than for the Tevatron (high pt). Looking at specific observables we 
see that, as already remarked in the case of Tm,g for all observables the symmetric scale 
variation uncertainties are considerably smaller than the full uncertainties, and we stress 
that only the latter are really indicative of the size of all kinds of neglected higher order 
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Figure 4: The normalised NLO+NLL matched distributions, for a range of event- 
shape observables at the Tevatron with pn > 200 GeV. For further details, see text. 
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Figure 5: The normalised NLO+NLL matched distribution, for a range of event-shape 
observables at the y/s = 14 TeV LHC with pn > 200 GeV. For further details, see text. 
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terms. 

Some final remarks concerns the NLO+NLL results for T^,??,, -Bt.tj and -Bvy,7^• As 
discussed in p3] and at the end of Sec. 13.1.21 recoil variables are more difficult to resum than 
other variables, because in the CAESAR approach, the NLL term g2{asL) has an unphysical 
divergence at small values of the observable. This difficulty is reflected in the substantially 
larger uncertainty bands for these observables than for the directly global variants and those 
with an exponentially-suppressed forward term. Among the recoil variables, the thrust 
minor and broadenings were the only ones for which an even partially acceptable result 
could be obtained. In order to obtain results for recoil variables of similar quality to those 
for the other observables requires a resummation of initial state emissions in appropriate 
Fourier transform variables, as done e.g. for the Drell-Yan pt resummation [75], mixed with 
a Sudakov type resummation, as was done for the DIS broadening [76]. This is beyond 
the scope of CAESAR. Another characteristic to be commented on is the spike for Bt^-r, 
and ByyTi near 0.37. We believe this could be related to a Sudakov shoulder [77] type 
phenomenon, and similar (though less pronounced) artifacts have also been observed in 
DIS event-shape distributions. Again, it is beyond the scope of CAESAR to resum the 
enhanced higher-order terms associated with these structures. 

4.5 Naive exponentiation of NLO 

In the previous Section we presented full NLO-I-NLL resummations for a range of event 
shapes. Both the NLO Monte Carlo calculation and the NLL resummation are highly 
CPU intensive and are usually both run across many CPUs. While the NLO part of 
the calculation is the most computer intensive, this is to some degree counterbalanced 
by the fact that many observables can be computed in the same NLO run. The NLL 
resummation on the other hand requires essentially a separate run with CAESAR for each 
observable. Altogether, a single combination of collider energy and event-selection cuts 
requires many CPU-years of computing time. NLO-I-NLL resummation also requires that 
the NLO total cross section and LO distributions be decomposed into flavour channels, 
and this information is not available in the public version of NLOJET++ (nor in most other 
public NLO codes). Furthermore, CAESAR is currently not public, the range of observables 
that can be resummed with CAESAR is not as broad as one might like (see Sec. I3.1.ip . 
and the matching procedure at hadron colliders is not as straightforward as in e^e~, as 
discussed in Sec. 13.21 For the above reasons, it is interesting to explore the possibility 
of obtaining predictions with accuracy close to NLO-I-NLL using publicly available NLO 
results only. For instance the following combination of LO and NLO integrated, flavour 
summed distributions, Si(t>) and S2(f), and the corresponding total cross sections, o"o and 
(Ti, as obtained directly from NLOJET++, 

^ + 5iM _ i pi(^)yj ^ (4.1) 

(To (To 2 \ ao / ' 

has the following properties: 



32 



0.1 




:«'kX I 1 


1 da 


)C/ ^^^^ 











NLL+NLO 

exp NLO ^rrm 



Tevatron, 1 .96 TeV 
Pti > 200 GeV, IVjetsI < 0.7, tic = ^ 



0.1 0.2 0.3 0.4 0.5 0.6 



Figure 6: Comparison of the NLO+NLL results with full uncertainties and exponentiation 
of NLO (as defined in eq. (14. ip ) with asymmetric scale uncertainties for the directly global 
thrust minor (top left), the total broadening with recoil term (top right) and for y^s 
(bottom left). Shown for the Tevatron with a 200 GeV cut on pti. 



• /(f) goes to 1 at f = fmax5 without O (a^) corrections; 

• the fixed order expansion of the corresponding differential distribution, up to relative 
order reproduces the normalized NLO differential distribution; 

• the formal accuracy is not even LL^: starting from order the terms a^L'^^ are 
only correct in the limit in which the event sample is dominated by a single colour 
channel. 

Though the method does guarantee any formal resummation accuracy, it is still instruc- 
tive to see how it fares in practice. We therefore show in fig. [Hla comparison of NLO+NLL 
matched results with full uncertainties and the naive exponentiation of NLO, as defined in 
eq. (14.11) with full scale uncertainties for a representative set of observables for the Tevatron 
with a 200-GeV jet pt cut (similar results hold also at the LHC and for other cuts as well 
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as other observables [S]). We see that for the directly global thrust minor (top, left) the 
exponentiation result is well-contained in the uncertainty band of the NLO+NLL matched 
result, suggesting that the naive exponentiation of NLO is indeed a quite reasonable pro- 
cedure (similar results hold in general for global observables). The same observation is 
true also for the total broadening with recoil term (top, right) but with one important 
difference. In this case the NLO+NLL uncertainty band at small values of the observable 
is divergent, signaling the breakdown of the resummation (as discussed in [23] and at the 
end of Sec. l3.L2p . The scale uncertainty band of the naive exponentiation does of course 
not account for this and is small for all values of the observable. Therefore for observables 
like the recoil event shapes, whose double logarithms do not fully exponentiate, the naive 
exponentiation of the NLO results can only be used as long as one is far from the divergence 
of the pure NLL resummation (whose position is one of the pieces of information provided 
by CAESAR). 

Finally, we show in fig. E] (bottom, left) how well the naive exponentiation does in the 
case of the three-jet resolution parameter with exponentially suppressed term. We see that 
in this case the naive exponentiation result is not contained in the full uncertainty band 
of the NLO+NLL resummation. This is true for the tail of the distribution, where there 
seems to be too little radiation suppression and, similarly, for the peak, whose position 
is slightly displaced to the left. This softer spectrum is a feature of other exponentially 
suppressed observables as well (though the effect is not as remarkable). 

Altogether it seems that naive exponentiation is a sensible procedure to extend the range 
of validity of pure NLO predictions. However, since it is not guaranteed to provide any 
formal logarithmic accuracy, one should rely on full NLO+NLL predictions for precision 
studies. In any case, we stress that before carrying out this exponentiation procedure, one 
should understand the basic soft /coUinear properties of the observable (be it with CAESAR 
or in whatever other way). 

4.6 Comparison with (matched) parton showers 

For most practical applications, it is far more convenient to use parton-shower Monte Carlo 
event generators, like Herwig [3] or Pythia [4\, or event generators merged with LO matrix 
elements, rather than a full NLL+NLO calculation. It can however be difficult to estimate 
the accuracy of these tools and the reliability of the error estimates that come with them. 
The purpose of this section is therefore to compare the NLL+NLO results with parton- 
shower based predictions (at parton level, in order to avoid non-perturbative corrections 
from hadronisation and the underlying event). 

We will start with Herwig (v6.5) events showered from exact tree-level matrix elements 
for 2 — >• 2, 2 — 7- 3 and 2 — )■ 4 partonic scatterings, as generated with Alpgen [78]. We use 
the MLM prescription [79] to avoid double counting between emissions generated in the 
hard 2 — )■ n scattering and those generated by the parton shower. This combination of 
parton-shower and tree-level matrix elements is the standard tool for many Tevatron and 
LHC predictions. 

To gauge uncertainties in the resulting matched samples we shall simultaneously vary 
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the renormalisation and factorisation scales in the tree-level matrix elements by a factor 
of two around their default settings in the MLM procedure (which are taken as in the 
CKKW procedure [80]). The MLM procedure also involves a separation scale between 
the region of phase-space to be accounted for by the tree-level matrix elements or by the 
parton shower. We take this separation scale to be 0.5pi,min when we look at event samples 
whose hardest jet has pt > Pt,mm and the angular distance for a jet and a parton to be 
matched is restricted to be AR < 1.05. The hard events have been generated with a pt 
threshold 0.4pf mm for all partons (constrained to have \y\ < 5), which must separated from 
each other by a distance AR > 0.7. In principle the MLM separation scale should also 
be varied in order to gauge uncertainties. However, the generation threshold should also 
be kept lower than the separation scale and the production of the 2 — )■ 4 tree-level event 
samples with the 0.4pt rnin threshold already turned out to have very low efficiency and 
would have become prohibitive with much lower a threshold. Therefore we will only show 
results with a fixed separation scale. The PDF that we used was CTEQ5L [81\, the default 
choice in Alpgen, but we also verified that the effect of switching to CTEQ6M (as used in 
our NLO+NLL calculations) was small. 

A comparison of the parton-level Alpgen+Herwig (Tree+PS) results with the NLL+NLO 
results is given in figs. [7] and [8|, for pp collisions at Tevatron and LHC, with a cut of 200 GeV 
on the transverse momentum of the hardest jet. The Alpgen-|-Herwig results are shown as 
red cross-hatched bands. The NLO-I-NLL results are shown as two bands: a blue hatched 
band whose width corresponds to the uncertainty from just the symmetric variation of 
renormalisation and factorisation scales; and a cyan, solid band corresponding to the full 
set of uncertainties represented in fig. [H 

Generally, there is reasonable agreement between the Tree+PS and the NLO+NLL 
results. One feature of note is that the Tree+PS uncertainty band is significantly narrow 
than the NLO+NLL band (even that with just symmetric scale variation). It is not imme- 
diately obvious that this truly reflects smaller uncertainties in the Tree+PS case, which, 
based as it is on LO calculations, would be expected to show larger uncertainties than 
the NLO+NLL prediction. We tend instead to interpret this as indicating that symmet- 
ric scale variation does not provide a good estimate of the true uncertainty on Tree+PS 
predictions]^ 

There does not seem to be a clear pattern to the cases where there are significant 
differences between the two kinds of predictions. For example, for the T±^g and the Tmin 
variables the Tree+PS predictions seem harder than the NLL+NLO results. Instead, for 
?/3_£, the Tree+PS results are generally softer. 

We also show in figs. M and [TU] a comparison between results obtained from different 
shower Monte Carlo event generators with and without matching, for the same subset of the 
observables at the Tevatron and the LHC with a minimum pt on the hardest jet of 200 GeV. 
Pythia 6.4 is shown both for the old (virtuality ordered) and new (transverse-momentum 
ordered) showers. All results are shown at parton-level, without multiple interactions (i.e. 



^■^It should be said that the uncertainty band on the Tree+PS prediction is considerably larger if one 
considers differential cross sections instead of normalised differential distributions. 
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Tevatron, 1 .96 TeV NLO+NLL (all uncert.) mtm 

Pti > 200 GeV, lyjgtsl < 0.7, r\Q = 1 NLO+NLL (sym. scale uncert.) K\\1 

PARTON LEVEL NO UE Alpgen + Herwig (partons) 

Figure 7: Comparison of the NLO+NLL results with matched Alpgen+Herwig results. 
The latter have just symmetric scale variation in their uncertainty bands, so we also in- 
clude subsidiary bands for NLL+NLO with just symmetric scale variation. Shown for the 
Tevatron with a 200 GeV cut on pn. 
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Figure 8: Comparison of the NLO+NLL results with matched Alpgen+Herwig results. The 
latter have just symmetric scale variation in their uncertainty bands, so we also include 
subsidiary bands for NLL+NLO with just symmetric scale variation. Shown for the ^/s = 
14 TeV LHC with a 200 GeV cut pa. 
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Tevatron, 1 .96 TeV NLO+NLL (all uncert.) mtm 

Pti > 200 GeV, lyjgtsl < 0.7, r\Q = 1 NLO+NLL (sym. scale uncert.) K\\1 

PARTON LEVEL NO UE Alpgen + Herwig (partons) 

Herwig 6.5 

Pythia 6.4 virtuality ordered shower (DW tune) 

Pythia 6.4 P( ordered shower (SOA tune) 

Figure 9: As is fig. [7] with, in addition, results from plain Herwig and Pythia showers. 
Pythia 6.4 is shown both for the old (virtuality ordered) and new (transverse- momentum 
ordered) showers. All results are shown at parton-level, without multiple interactions (i.e. 
no underlying event). Shown for the Tevatron with a 200 GeV cut on pn. 
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LHC, 1 4 TeV NLO+NLL (all uncert.) 

Pti > 200 GeV, lyjgtsl < 1 , tiq = 1 .5 NLO+NLL (sym. scale uncert.) 

PARTON LEVEL NO UE Alpgen + Herwig (partons) 

Herwig 6.5 

Pythia 6.4 virtuality ordered shower (DW tune) 

Pythia 6.4 P( ordered shower (SOA tune) 

Figure 10: As is fig. |8] with, in addition, results from plain Herwig and Pythia showers. 
Pythia 6.4 is shown both for the old (virtuality ordered) and new (transverse- momentum 
ordered) showers. All results are shown at parton-level, without multiple interactions (i.e. 
no underlying event). Shown for the y/s = 14 TeV LHC with a 200 GeV cut on pti. 
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Herwig 6.510 

Pythia 6.412 DW - - - 

Pythia 6.421 PerugiaO 

Pythia 6.421 SOA 

Pythia 6.421 Pro-pTO 

Tevatron, 1 .96 TeV 
> 200 GeV, |yj3,3| < 0.7, = 1 



Figure 11: The differential distribution for y-^ g at the Tevatron, pn > 200 GeV, computed 
at parton level with Herwig and Pythia with both the old shower (DW) and the new shower 
with different tunes. 

no underlying event). For reference we show also the result of NLO+NLL resummation 
with symmetric scale variation uncertainties (the use of the full band would complicate the 
plots excessively). 

In general Herwig's angular-ordered shower and Pythia's virtuality ordered (old) shower 
give results that are quite similar (or slightly harder) to the full matched results, with 
deviations visible in some cases, e.g. for at the LHC. 

It is perhaps surprising that unmatched, plain parton shower results, are often harder 
and sometimes even closer to the NLL+NLO band than the Tree+PS matched ones, this 
is particularly evident for y^^g at the LHC. This is an unexpected result as the motivation 
for carrying out Tree+PS merging is that parton showers are unable to reproduce the 
structure of hard large-angle emissions. 

What is also evident from figs. M and dO] is that there are big discrepancies between 
the newer, transverse-momentum ordered shower in Pythia 6.4 (in the SOA tune) and, 
essentially, everything else. These distributions appear to be significantly softer than those 
from other parton showers, with the difference most visible in the case of the variables, 
and inconsistent with the NLL+NLO calculation. 

It is therefore useful to try to understand whether the origin of the discrepancies lies in 
the new shower or in the tuning of the Pythia parameters. To further probe this issue, we 
show in Figure [11] a comparison among plain Herwig, virtuality ordered (DW) Pythia and 
different tunings of the new transverse- momentum ordered shower, SOA [82], as used above, 
and two more recent tunes, PerugiaO |83j Pro-ptO [SI], (shown with version 6.421; version 
6.412 yields identical results). Of these the two more recent pj-ordered tunes, Pro-ptO gives 
results very similar to SOA, while PerugiaO is closer (though not identical) to the Herwig 
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and virtuality-ordered Pythia results. The conclusion to be drawn from these results is 
that for transverse-momentum ordered showers, the shower parameters can have major 
implications for the reliability of the results and a consensus has yet to emerge among 
current tunes for the choices of these parameters. 

5 Non-perturbative effects 

So far we discussed only perturbative effects, however, before any comparison to data can 
be done, non-perturbative effects have to be included. As in e^e~ annihilation, there are 
non-perturbative effects due to hadronisation, i.e. related to the transition of partons to 
hadrons. In hadron-hadron collisions there are also interactions between the two beam- 
remnants, the so-called underlying event (UE). Both effects are suppressed by inverse 
powers of the hard scale pt of the high-pt scattering. For event shapes the dependence 
is linear in 1/pt |H5], while for jet-resolution parameters the dependence is even more 
suppressed, as will become evident also from the plots presented in this section^ 

At hadron colliders, analytical predictions for non-perturbative (NP) effects are avail- 
able only for a limited number of jet-observbales [851 EHl [871 [HH]. A more general way to 
estimate those effects is to use event generators, such as Herwig [3]or Pythia [4], which 
can be run at parton- or hadron-level with or without underlying event. The default 
Pythia underlying event model includes multi-parton interactions, while Herwig needs to 
be interfaced to Jimmy [89] to have a realistic modelling of the underlying event. 

In Figs. [T2lfT5l we compare parton level and hadron level results without and with 
UE for the set of observables discussed previously for the low- and high-p^ samples at the 
Tevatron and at the LHC, as obtained with virtuality ordered Pythia 6.4 (DW tune). 

As far as hadronisation corrections are concerned, one notices immediately that for 
event shapes these effects are quite large at the Tevatron for the pti > 50GeV sample. 
They systematically shift the distributions to the right and distort them (mostly squeeze 
them). As expected these effects decrease considerably at pa > 200GeV. Going from 
the Tevatron to the LHC, keeping the pti > 200GeV cut, hadronisation correction are 
comparable, while again they decrease when going to the pn > lOOOGeV sample, where 
they are completely negligible. Since the majority of events in a sample will have jets 
with Pt close to the Pt-cut, this patten confirms the expected 1/pt scaling of hadronisation 
corrections. For distributions hadronisation effects follow the same pattern but are much 
smaller, and are already very small at the Tevatron for pn > 50GeV. 

The effect of the underlying event on these distributions is quite different. For event 
shape distributions at the Tevatron for pn > 50 GeV, the UE broadens significantly the 
distributions and moves them systematically to the right. For fixed center of mass (c.o.m.) 

A further potential non-perturbative effect is that due to in-time pileup, the additional (usually) soft 
pp collisions that take place during the same beam crossing as the hard pp collision of interest. Its impact 
can largely be eliminated by considering only events with a single primary interaction in the beam crossing. 
Given the huge cross-sections for the event selections outlined in table [21 this should not pose too much of 
a problem except, possibly, at the higher LHC pr cut. 
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Figure 12: Comparison of parton- level, hadron- level without UE and hadron level with 
UE, for selected event-shape distributions, as obtained with Pythia 6.4 (DW tune). Shown 
for the Tevatron with a 50 GeV cut on pti- 
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Figure 13: As in fig. [121 but for the Tevatron with a 200 GeV cut on pti 
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Figure 14: As in fig. HH but for tlie = 14 TeV LHC witli a 200 GeV cut on pa. 
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Figure 15: As in fig. [121 but for tlie y/s = 14 TeV LHC witli a 1 TeV cut on pn. 
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Tevatron, 1 .96 TeV Herwig + Jimmy (ATLAS tune) 

Pti > 200 GeV, lyjgtsl < 0.7, ^Q = 1 Pythia DW 

HADRON LEVEL WITH UE Py*'^'^ SOA 



Figure 16: Hadron-level results, including underlying event, for selected event-shape dis- 
tributions, as obtained with Herwig 6.5 -|- Jimmy 4.31 (ATLAS tune, as given in text) and 
Pythia 6.4 with 2 tunes, DW (DWT would be identical) and SO A. Shown for the Tevatron 
with a 200 GeV cut on pa. 
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Figure 17: Hadron- level results, including underlying event, for selected event-shape dis- 
tributions, as obtained with Herwig 6.5 + Jimmy 4.31 (ATLAS tune, as given in text) and 
Pythia 6.4 with 3 tunes, DW, DWT and SOA. Shown for the = 14 TeV LHC with a 
200 GeV cut on pa. 



energy, the UE decreases with pt, but contrary to the hadronisation corrections, increasing 
the c.o.m. energy, at fixed pt results in an increased UE activity. As for hadronisation 
corrections, UE effects are much smaller for distributions than for event shapes. This 
fact means that at sufficiently high pt one can compare perturbative predictions directly 
to data, without additional NP corrections. This also makes 7/3 distributions (in particular 
the global version) suitable for direct tuning of shower parameters. We note also that 
different event shapes have different NP sensitivities, broadenings seem to have smaller 
corrections, while masses and thrust distributions tend to have larger ones. Therefore the 
latter seem better suited to study NP effects and to tune models of hadronisation and 
underlying event. 

To address this last issue further, we show in Fig. [16] how different Monte Carlo showers 
and tunes to the same Tevatron data differ from each other for the same set of observables. 
Specifically we use Herwig+JimmyEl Pythia's virtuality ordered shower with the DW tune 
and Pythia with the p^-ordered shower (SOA tune). It is noticeable how the DW, SOA and 
Herwig+ Jimmy tunes differ (sometimes substantially), despite the fact that all have been 
tuned to Tevatron data. For SOA in particular this is not really surprising: if perturbative 
predictions already differ substantially, so will full results at hadron level. However it 
is nevertheless instructive, because it illustrates to what extent hadron-level event-shape 
distributions can help constrain perturbative aspects of the shower. 

Finally, in Fig. [T7]we show what happens for the same pt cut at the LHC Discrepancies 
between Herwig and Pythia survive (but are maybe reduced). In addition to the DW 
Pythia tune, we also show DWT (which was identical at Tevatron energies) and see sizable 
differences between them. All this suggests that event shapes have significant scope for 
tunes of event generators. 



6 Multi-jet limit 

One common use of event shapes is to distinguish two different classes of multi-jet events: 
those in which the jets cover phase space quite uniformly, as in multi-body heavy-particle 
decays; and those in which the jets are relatively collimated in few bunches, as induced 
by the coUinear singularities of massless QCD multi-particle emission. For example, the 
ATLAS [51] and CMS [52j discussions of prospective analyses mention the use of event 
shapes, most notably the transverse sphericity (see below), in physics studies that range 
from tt analyses to searches for supersymmetric particle decays and black-hole decays (yet 
other applications include hidden- valley searches [nil[S2])- 

The purpose of this section is to compare various event shapes' ability to distinguish 
characteristically different event topologies. Firstly we shall examine to what extent current 
event shapes are able to distinguish symmetric 3-jet topologies from symmetric multijet 
topologies. The main finding will be that they discriminate principally between 2-jet pencil- 
like events and symmetric events, regardless of the number of jets in the latter. We shall 

15 The parameters used are PRSOF^O, PTJIM= 2.8(s/1800 GeVf -^^'^ GeV, an inverse (anti)proton 
radius of JMRAD(73)=JMRAD(91)= 1.8 GeV and CTEQ6L1 68J PDFs, as per the ATLAS tune in f90]. 
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then study the robustness of the identification of symmetric events: both with respect to 
parton showering and to the orientation of the multijet system. The results from these two 
studies will then lead us to propose event shapes that should have enhanced sensitivity to 
the symmetric multijet limito 



6.1 The transverse sphericity 

One event shape that we have not considered so far is the sphericity. Since it is by far the 
most widely used for discriminating symmetric multi-jet topologies, let us briefly examine 
its properties. It is defined in terms of the eigenvalues Ai > A2 of the transverse-momentum 
tensor: 

M _ I PxiPyi \ qpheri _ ^-^2 (a ^\ 

^^y-z^[ p^^Py, pI ) ' ^^.^ = aTTa; • ^^■^> 

It has the property that it tends to 1 for events with circular symmetry in the transverse 
plane, and is for pencil-like events. However the appearance of a sum of squared momen- 
tum components in M^y makes this observable collinear unsafe, as is the case [M] for the 
related variable in e^e~: for example, if a hard momentum along the x direction is split 
into two equal collinear momenta, then their combined contribution to ^ ■ p^.^ will be half 
that of the original momentum. Therefore collinear splittings change the sphericity by a 
factor of order 1. One consequence of this is that it is impossible to make perturbative 
predictions for the sphericity beyond leading order. Another consequence is that parton 
showering and hadronisation significantly alter the value of the observable, limiting its 
ability to discriminate different topologies (at least when the input momenta are particles; 
often it is jets that are used as inputs). We shall see this exphcitly in section | 



6.2 The circular Umit 

The simplest instructive study that comes to mind for event shapes intended to distin- 
guish symmetrical multi-jet events from dijet events is to examine their values V^(A^) 
for perfectly symmetrical planar transverse events with varying numbers of momenta. 
Pi = %^(cos ^, sin ^, 0, 1) for z = 1 . . . N. This is illustrated in fig. CHI For uniformity 
the results have been normalised to the value Vdrc obtained for perfectly circular planar 
events (A^ — )■ oo), as given in table [31 

The only two observables with a monotonic (and trivial) behaviour are g and Fg. 
The remaining ones have been grouped into the left and right hand plots according to 
whether they peak for n = 3 (thrust-like) or n = 4 (broadening-like) and one sees that 
the perfectly circular limit does not give the largest value for all observables. Perhaps 
the most interesting observation from these plots is the modest difference between the 

It is worth noting also the event-shape type observable proposed for BSM searches [93], whose aim is 
not to distinguish multijet "hedgehog" topologies from dijet topologies, but rather to be sensitive to the 
hadronic structure of BSM signal events with large transverse missing energy (e.g. i?-parity conserving 
SUSY) , but without making explicit use of the measurement of missing transverse energy. 
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Figure 18: Values of various observables for events with N momenta arranged symmetri- 
cally in the transverse plane. 
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Table 3: Kirc? the values of various observables in the transverse circularly symmetric limit. 
The events have been chosen to be planar — the variables other than ps^c and Bt,c are 
however insensitive to the longitudinal components of the momenta. 



3-particle and fully circular events — thus they are sensitive to the absence of a unique 
preferred transverse direction, but not to the overall degree of symmetry of the event. 



6.3 Collinear safety and showered events 

Let us now ask the question of how much the collinear unsafety of S^^" matters in prac- 
tice. To investigate this we have taken a number of 2 — )■ 3 partonic events and showered 
them with Herwig, using the "inclusive" MLM prescription [78] to reject events in which 
the showering introduces a fourth, harder jet, or other strong modifications of the event. 
Figure [19] shows the distribution of S^^", Fg and Bt^c for two kinds of input 2 — )• 3 
partonic event, which are both planar with all particles at rapidity y = 0: 





Event 1 ( 


generic) 


Event 2 (Mercedes) 


Ptl 


= 828 GeV, 


01 = 


Ptl 


= 666 GeV, 0i = 


Pt2 


= 588 GeV, 


02 = 37r/4 


Pt2 


= 666 GeV, 02 = 27r/3 


Pt3 


= 588 GeV, 


03 = -37r/4 


Pt3 


= 666 GeV, 03 = -27r/3 



One clearly sees that the collinear unsafe S^^" has much less discriminating power between 
the two events than do Fg or Bj-c (or for that matter any of the other event shapes that 
were shown in fig. [T5]) . Furthermore it is clear for the Mercedes event that the peak at 
gphen _ Q Yia,s little connection with the expected 3-parton Mercedes value of 1. In 
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Figure 19: Distribution of event shape values after showering and hadronisation for the 
"generic" and Mercedes input 3-parton events (for further details of the event generation, 
see text). The arrows indicate the values for the 3-parton events. Small overlap between the 
two distributions and good correspondence with the arrows are signs of a good observable. 



contrast, the distributions for the other two observables are peaked close to the expected 
values (indicated by the arrows). This should of course be of no surprise given the collinear 
unsafety of S^^^". However, in view of the latter's widespread current use (albeit with jets, 
rather than particles, as inputs), we feel that the point is worth noting. 

6.4 Impact of event orientation 

One of the interests of event-shape studies is in identifying massive particle decays. Most 
of the event shapes above have the counterproductive characteristic that they give very 
different results for particles that decay with just transverse components (in the particle's 
centre of mass) or with both longitudinal and transverse components. To illustrate this, 
we take the generic event given above and rotate it by 11/2 around the axis of particle 1, 
giving pt2 = Pt3 ~ 416 GeV, 02 = 03 = tt and rapidities 1/2 = — 2/3 — 0.88. We shower it, 
as explained above, and the resulting distributions for three event shapes (normalised to 
their values in the circular limit) are shown in fig. [201 

For S^g" and Fg there is a large difference between the distributions for the generic and 
rotated-generic events (and similarly for e.g. 8^^^^° , T±^g and T^in,g)- For the broadening 
in contrast, which we recall involves both the y and dispersions of particles with respect 
to axes in each of the two central half-regions, the generic and rotated-generic events give 
rather similar distributions. A similar phenomenon occurs with the invariant masses of 
those regions, in that masses too are sensitive to both directions of dispersion, though their 
intrinsic rotational invariance is in part spoiled when one normalises to Q±fi as in eq. fl2.17p . 
The rotational invariance is probably in part the origin of the usefulness of "cluster-masses" 
in the context of hidden- valley studies [92j . Note however that masses are significantly more 
sensitive to (initial-state) forward semi-hard radiation than are broadenings. 
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Figure 20: Distribution of event shape values after showering and hadronisation for the 
"generic" and rotated generic input 3-parton events, as described in the text. Observables 
that have similar distributions for the two sets of events are likely to be more effective for 
identifying massive-object decays. 



6.5 Increasing sensitivity to the spherical hmit 

As is clear from figure [181 none of the event shapes above are particularly effective at distin- 
guishing truly spherical events from simpler multi-jet topologies, like symmetric transverse- 
planar events. 

What one has in mind when discussing spherical events is that they have significant 
"volume", symmetrically distributed around the event. One way of obtaining sensitivity 
to this is to consider the following matrix, separately for the up and down central regions 
of an event: 

2fcC(7 

where Aynj = t/i — yc^ and A(j)iu = (pi — (pcu, and similarly for the central down region, 
Cd- The eigenvalues Xui > \u2 of Mu have the property that Xui is non-zero if there 
are two non-coUinear particles in the hemisphere, while Xu2 is non-zero if there are three 
non-coplanar particles in the hemisphere. The observable 

= min(A(72, Ad2) , (6.3) 

which we name "supersphero" , is therefore non-zero only if there are 3 non-coplanar par- 
ticles in each of the hemispheres of the event — i.e. for events that truly bear some 
resemblance to spherical events. 

For a perfectly spherical event the two eigenvalues in each hemisphere are Ai = 7r^/24 ^ 
0.411 and A2 = tt^/S - 1 ~ 0.2340 

The 5*6 observable, in terms of its use of eigenvalues of a 2 x 2 matrix, relates of course 
to the Fg-parameter of eq. (12. 4p and to the event shapes studied for boosted top-quark 

^"^ One may balance the two eigenvalues more closely for example by replacing Ayiu — > Aymil + Ayfu/G) 
(i.e. with the first two terms of the expansion of smhAynj) in eq. (|6.2I) . though there is a significant degree 
of arbitrariness in this choice. 
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identification [261 12H]- The latter's use of a matrix defined in tlie plane transverse to a jet 
is actually quite similar to our use of a matrix defined in a central half-region. 

A detailed study of the 5*6 observable would benefit from comparisons of high jet- 
multiplicity QCD samples and multijet samples from new-physics scenarios. Such a study 
is beyond the scope of this paper, but would, we believe, be of interest. 

7 Summary of main results 

Given the length of the paper, and the fact that we have addressed quite a range of issues, 
we find it useful, before concluding, to summarise here the main results of the paper. 

There are a number of reason why event shapes provide a powerful laboratory for 
studying of a range of aspects of strong-interaction physics at hadron colliders. From an 
experimental point of view, cross sections for the QCD (dijet) events on which one car- 
ries out event-shape studies are very large both at the Tevatron and at the LHC This 
means that high-statistics event samples are already available at the Tevatron and can 
be expected early on at the LHC. Since event shapes are defined as dimensionless ratios 
of combinations of hadron-momenta, and since their differential distributions are also di- 
mensionless, many experimental uncertainties are reduced. From a theoretical point of 
view, one of the attractive characteristics of event-shape studies is that different variables 
provide complementary sensitivities to a broad variety of features of hadronic events, such 
as the topology of the final state, the nature of initial and final-state jet-fragmentation, 
hadronisation, and the underlying event. 

The above points motivated us to study many hadron-collider event shapes in a sin- 
gle context. We exploited the automated NLL resummation procedure implemented in 
CAESAR to obtain next-to-leading logarithmic resummed distributions matched to next- 
to-leading order exact predictions from NLOJET++ for a large number of event shapes at 
the Tevatron and at the LHC. 

The matching procedure is conceptually simple, but technically involved, as discussed 
in detail in Sec. 13.21 One issue is that the best logarithmic accuracy achievable once NLO 
and NLL predictions are available, namely in the expansion of the integrated 

distribution, NNLL^, can be obtained only if the NLO code provides full information 
about the flavour of all incoming and outgoing partons. This decomposition into flavour 
channels is not present in the publicly available version of NLOJET++, so we used the 
extended version developed in [94] in order to extract this information. We also needed 
to use the fiavour-fct jet-algorithm of ref. [S5 to map the flavour of 2 — )■ 3 events into 
that of an underlying 2 — )■ 2 Born-like event. Additionally we needed the order-by-order 
expansion of PDF evolution, which was provided by hoppet [61]. The computing effort 
should also not be neglected: our directory of resummed results contains O (10000) files, 
and we estimate that several tens of years of CPU time have gone into the NLO and 
NLL calculations used here. Part of this complexity stems from our choice to consider 
several different classes of uncertainties associated with uncalculated higher-order terms: 
those from separate variation of renormalisation and factorisation scales; redefinition of 
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the argument of the logarithm being resummed (X-scale); and two choices of schemes for 
combining (matching) the NLO and NLL results. 

Among the questions we asked was whether this considerable complexity is needed. It 
turned out that the flavour decomposition had only a modest effect (cf. Appendix Rl). 
We also found that a simple exponentiation of the NLO result, as presented in Sec. 14.51 
not even correct to LL or LL^ accuracy, comes remarkably close to reproducing most of 
the NLO+NLL distributions (albeit not close enough that one would forgo NLO+NLL if 
it is available). One interpretation is that the large amount of radiation that comes from 
the 4 Born legs in a 2 — )■ 2 process causes event-shape distributions to be dominated by 
regions where the logarithm that is being resummed is not all that large. Note, however, 
that plain (unexponentiated) NLO predictions are very inadequate substitutes for the full 
NLL+NLO result and their uncertainty bands are misleadingly small. 

We studied three generic classes of event shapes: the directly global ones, those with 
exponentially suppressed terms and those with a recoil term. The definition of the observ- 
ables is recaleld in sections. 12. 12. 2[ and 12.31 respectively. While stable numerical results 
could be obtained for observables belonging to the first two classes, for the last of these 
it was sometimes impossible to obtain numerically sensible results. This is in part due to 
cancellations among contributions from multiple emissions in the recoil term, which cause 
the resummation provided by CAESAR to have a divergence at small w , as explained in 
Sec. 13.1.21 It is also due to structures in the middle of the physical region, akin to Su- 
dakov shoulders [77], which would require an additional resummation. Such shoulders are 
visible e.g. for the broadenings with recoil term in Figs. H] and |5l These observables are 
also challenging experimentally because the measurement of the recoil term is affected by 
cancellations between large transverse momenta of the two hard jets. 

A question mark that hangs over NLL resummations is that of coherence-violating 
logarithms (CVL, referred to as super-leading logarithms (SLLs) in the context of interjet 
energy flow) [39lll0], terms potentially starting at afL^, related to a violation of coherence, 
whose validity was a crucial assumption in the resummations of [31]. There is a risk 
that this could therefore invalidate our claim of NLL accuracy for some observables. We 
investigated this point in Sec. 13. 3^ and found that the answer depends critically on the 
ordering parameter used in the calculation of the SLL terms. If, as in [3S1 SO] , one makes 
the assumption that the ordering parameter is transverse momentum, then the claim of 
NLL accuracy breaks down for our "exponentially suppressed" class of observable (not 
for the others), while NNLLs remains valid for all observables. If one instead assumes 
virtuality ordering, then both NLL and NNLL^ accuracies should be valid for all our 
observables. This highlights the importance of understanding the question of ordering 
for SLLs, which also affects the coefficient of the a^L^ terms in [391 HQ] and probably 
requires that one go beyond the eikonal approximation that was used there. Nevertheless, 
practically we tend to believe that SLLs will not seriously affect our results, one reason 
being that we still retain NNLL^ accuracy. 

Turning to our phenomenological results, a feature common to all observables is that the 
shape of the distributions is strongly influenced by the ratio of quarks to gluons among the 
incoming partons. This is because the double-logarithmic Sudakov exponent, responsible 
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for the position and width of the peak of the distribution for each underlying subprocess, 
is determined by the total color charge of the hard emitting partons. Event samples 
dominated by gluon scattering (Tevatron with pa > 50 GeV, LHC with pu > 200 GeV) 
have broader distributions than those dominated by quark scattering (Tevatron with pa > 
200 GeV, LHC with pa > 1000 GeV). This is evident e.g. in Fig. E] in the case of our 
representative observable ^ and is discussed in Sec 14.31 We remark that dijet event- 
shapes at hadron colliders are the first case in which a change in a kinematical cut modifies 
the double logarithmic behaviour of the event-shape distributions. This would not be the 
case for event shapes in hadron-coUider processes such as Drell-Yan production, or Z+jet 
or H^-Fjet. 

In the absence of data on the event shapes discussed here, one of the interesting uses 
of our NLL+NLO results is to compare them to the results of two Monte Carlo parton 
shower programs, this is discussed in Sec. 14.61 We considered Pythia 6.4 and Herwig 6.5 
both without and (in the case of Herwig) with matching to multi-parton tree-level matrix 
elements (Alpgen, MLM prescription). The quality of the agreement between plain parton 
showers and the resummations depends significantly on the quark/gluon admixture: in 
quark-dominated event samples it is often adequate, while in gluon-dominated samples 
it is somewhat poorer. This may be a reflection of the extensive tuning of quark parton 
showers carried out with LEP data, while gluon parton showers have seen fewer constraints. 
The importance of tuning parton showers in a context with incoming beams is highlighted 
particularly strongly by the results of the newer pt-ordered shower in Pythia 6.4. In two 
tunes, SOA and Pro-PtO, the agreement both with NLO+NLL and with other showers is 
quite poor; in the PerugiaO tune it improves, as can be seen from Fig. [TTl 

One might expect that supplementing parton showers with matching to multi-parton 
tree-level events (Tree+PS) should improve the agreement with NLO+NLL results. This is 
the case only for some of the observables. We also examined the impact of (simultaneous) 
renormalisation and factorisation scale variation on the Tree+PS results and found that 
it leads to an uncertainty estimate that is far smaller than the actual differences between 
Tree+PS and NLO+NLL results, as can be seen in Figs. [7] and [Hll^ This should not be 
surprising: in the NLO+NLL calculations simultaneous scale variation represented only a 
small part of the full uncertainties. Questions that remain open therefore are whether in the 
Tree+PS approach uncertainties can be more faithfully estimated if one examines further 
"handles" (independent scale variation, matching scale, etc.), and whether we would have 
reached similar conclusions with other matching schemes (e.g. CKKW) and programs. 

From a non-perturbative point of view, we estimated both hadronisation and UE correc- 
tions using Monte Carlo event generators, as discussed in Sec. O As expected, hadronisa- 
tion corrections decrease when increasing the pti-cut on the jets. They are fairly negligible 
with cuts of the order of 200 GeV both at the Tevatron and at the LHC, as can be seen in 
Figs. [13] and [HI For lower pt cuts, they shift the distributions to the right and, for some 
observables they squeeze them, see e.g. Figs. [121 For jet resolution parameters {y^^g and 

^®More precisely: it was significant on the 2, 3 and 4-jet differential cross sections, but mostly cancelled 
in the normalised event-shape distributions. 
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y-i^g) hadronisation effects are always small, just a few percent correction for pti > 50 GeV 
at the Tevatron, much smaller in all other cases. These observations are consistent with 
the experience obtained from e~^e~ and DIS event-shape studies. 

As concerns the UE, there are observables for which it has a sizable effect even at 
Pti > 1 TeV, most notably for the thrusts and jet-masses, as can be seen in Fig. [151 This 
means that these event shapes are particularly good for tunes of the UE. Jet-resolution 
parameters are the only observables for which the UE effects remain consistently small (a 
few percent for the lower p^-cut samples, even smaller for the large pt ones). They are 
therefore well suited for tunes of perturbative parameters of showers and in general for 
perturbative studies. 

Finally, in Sec. |6] we examined how well event shapes can discriminate QCD-like two- 
jet events from BSM-like multi-jet events, and how robust this discriminating power is 
with respect to parton shower (radiative) corrections. In general we find that event-shapes 
discriminate well between events with two or more than two jets, but they do not discrim- 
inate well between three or any large number of jets: the value of event-shapes does not 
even increase monotonically with the number of jets for symmetric events, see Fig. [18] in 
Sec. 16.21 On the other hand it is possible to design new event shapes, which start with 
six jets in the final state, as is the case for our "supersphero" 5*6 event shape defined in 
Sec. 16.51 We believe these might be particularly promising for extracting new-physics sig- 
nals that involve relatively isotropic events with high jet multiplicity. Other considerations 
that we examined in Sec. 16.31 include how well event shapes retain their discriminatory 
power after parton showering (the collinear-unsafe, but widely used transverse sphericity, 
whose definition is recalled in Sec. 16. H is particularly poor in this respect); and also their 
sensitivity not just to transverse event structure, but also to longitudinal event structure 
(the broadenings do well at treating both on an equal footing). 

8 Conclusions 

In this article we have shown the first NLO+NLL (NNLL^) predictions, with full uncer- 
tainty bands, for hadronic observables at pp and pp colliders. We opted to make these 
predictions for event shapes in the context of dijet production, bringing together calcu- 
lations with CAESAR and a specially adapted version of NLOJET++, despite the fact that 
the NLO-I-NLL matching is technically more challenging than for event shapes in other 
hadron-collider processes such as Drell-Yan [24] or W/Z+jet [22] production. 

Several properties of the dijet process motivated our choice: it involves both initial 
and final-state partons; it offers the freedom to vary the proportion of quarks and gluons 
involved in the Born process, through the cut on the hard jets; when that cut is placed 
at moderate pt, dijet production involves a substantial gg — )■ gg scattering component, 
offering the most accessible example of a gluon-dominated process; and the cross sections 
imply large event samples. 

Comparisons of our results with parton-shower Monte Carlo predictions revealed ade- 
quate agreement for historic showers (Herwig 6.5, virtuality-ordered Pythia 6.4) in quark- 
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dominated cases, while the showers were generally too hard in gluon-dominated processes. 
Some common tunes of the newer, p^-ordered shower in Pythia 6.4 fared noticeably worse 
than the historic showers. We also examined one framework for matching to multi-parton 
tree-level matrix elements (MLM matching of Alpgen+Herwig 6.5). Though it led to some 
improvements, it was not immediately sufficient to bring about systematic agreement with 
the NLO+NLL results. These findings illustrate how event shapes can provide substantial 
input to the quest of understanding perturbativc QCD at hadron colliders. 

At hadron level, some event shapes are subject to significant non-perturbative correc- 
tions from hadronisation and the underlying event. We saw this to be the case, for exam- 
ple, for the thrusts and jet masses, while other observables, notably the ys variants, were 
largely unaffected by non-perturbative effects. Studying a broad range of event shapes, as 
done here, therefore provides complementary information on QCD phenomena at hadron 
colliders at many different physical scales. 

Event shapes are of interest not just for constraining QCD dynamics, but also for 
discriminating BSM-like multi-jet topologies from more QCD-like events. There are many 
interesting questions to ask about event shapes in this context. Some that we addressed 
here include their robustness to parton showering (the widely used transverse sphericity 
fares poorly), their sensitivity to longitudinal versus transverse event structure and their 
behaviour in the high jet multiplicity limit, where new dedicated event shapes, like the 
supersphero variable introduced here, can have particular advantages. 

These first steps of ours in exploring the phenomenology of event shapes at hadron 
colliders open a window onto a broad range of possible new studies, both theoretical and 
experimental. We look forward to their future development. 
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A Cross-checking fixed order and resummation 

Part of the value of having separate resummed and fixed-order calculations for event-shape 
distributions is that they provide cross-checks as to the validity of each of the approaches. 
This check is usually performed by a comparison of the exact fixed order results Si(v) in 
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eq. (13. 5p with the expansion of the resummed result T,r{v) from section 13.11 At small v 
the two results should differ order-by-order only by terms suppressed by powers of v or 
by logarithmically enhanced terms that are neglected within the resummation accuracy of 

At order as, the distribution S(t>) of eq. (13. 2p has the expansion at small v 



a 



(A.i: 



Sea 



where -ffim is the coefficient of L"^, has the dimension of a cross section and implicitly 
contains a" (notice that fflm is of order a""^^). A NLL resummation predicts ifjg'* ^^'^ 
h['^\ while hI^^ is obtained from the coefficient (agCj^^) of eq. (I3.18P by summing over all 
subprocesses S corresponding to the same colour channel a, as indicated in eq. flA.ip . The 
constant i/Jo^ can be extracted from the exact fixed cross sections cr^"^ and S^"''(f), defined 
as in eq. (13. 6p . as follows: 

Hil> = aS'^) + hm - ( Ht2L' + H^^ L)] . (A.2) 

v^Q L V /J 

Fig. [2Th shows the prediction for the differential distribution v[d'Li^j.{v) / dv\ obtained from 
eq. flA.ip . compared to the exact result v[(nii{v) / dv\ from NLOJET++, for the total trans- 
verse thrust r_L The two distributions agree at small v. 

Since fig. [2Tk contains large logarithms, a better visual constraint can be obtained by 
plotting the difference between Si(f) and its logarithmically-enhanced part H12LP' + HuL, 
which should go to a constant at small v, and indeed does. By performing this exercise 
separately for each colour channel one can obtain the h'^ individually, and can also verify 
that T!f^^^'^\v) vanishes for small v. From h''^ one can extract the colour-decomposed 
average coefficient constant {asC'f'^) = h[q^ /cXq"'\ 

The coefficients {agCl""^) obtained in this way are not precisely the ones that multiply 
the resummed distribution according to either of the two matching procedures described in 
section 13. 2[ because there one resums not logarithms of v but of a rescaled quantity Xyv, 
eq. (I3.24p . To get an idea of the size of the O (a^) term as it is relevant in the matched 
resummations, instead of plotting J2i{v) — {Hi2L'^ + HuL) , in fig. fUb we plot the difference 
between Si(f) and the distribution 

i:r,iiv) = H^2L^ + HuL = H^2L^ + HuL + H[o, L = \nJ—, (A.3) 

where Xy is the constant Xy of eq. (I3.25P computed for the reference Born configuration 
used for the analysis of the event-shape properties in CAESAR (two hard jets in the centre-of- 
mass frame with an angle 6* with respect to the beam corresponding to cos^^* = 0.2). The 
constants Him and H[q are defined in terms of the Him so as to give equality between the 
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middle and right-hand sides of eq. (]A.3|) . One observes that the difference ) — T,r,i{v) 
in fig. [2Tb (normahsed to ctq) goes, as expected, to a constant. That constant should be 
of order a^, whereas numerically it is O (1). However, we also know from table [2] that the 
order corrections can come with large coefficients. 

Given the {asC["'^) and the corresponding NLL resummations /'•"^(f), one can predict 
the NNLLs terms in the a<j expansion of i.e. terms a^L^ with n — 2 < p < n. 

Specifically, to second order in a^, we have 

T.2{v) = H2aL* + H2zL' + H22L^ + + (A.4) 

and to NNLLs accuracy we should control H2i, H23 and H22- To see that this is the case 
we compare E2(w) to the resummation prediction for the modified integrated distribution 

^vA^) = H2aL^ + ^23^' + H22L^ = H2aL^ + H2zL' + H22L' + H'^^L + H'^^ , (A.5) 

where again the constants H2m and ifgm defined so as to given agreement between the 
middle and right-hand sides of eq. ( 1A.5I) . Fig. [2Tb shows the exact second-order differential 
distribution f[(iS2(f) /(if] compared to t'[c/S2,.r('y)/'^'y] obtained from eq. ( lA.Sp . Again one 



sees good agreement, which is more readily verified by examining the difference between 
the two distributions, fig. [2Tl d. which is supposed to be (and is) flat (the constant results 
from differentiation of the H21L term in eq. (IA.5I) ). We also include the result that is 
obtained (lower points with errorbars) if one does not carry out the colour decomposition 
for {agCf^), but just computes (asCi) = Hiq/ctq. This gives rise to a different expansion, 
S^2('^)5 whose coefficient of is different from that of S,.^2(f)- For r^^g one notices that 
the corresponding difference between the exact result f [(iS2(f )/(if] and the distribution 
f [(iEg r(f )/(if] exhibits a hint of a slope at small r^^, indicating a missing alL"^ term in 
E^,,(f). 

Fig. [22] shows the same comparison of fig. [21] for the global thrust-minor Tm,g- In this 
case one is not able, within errorbars, to see any difference between a resummed prediction 
containing {usCf^), giving the correct H22-, and one based on (a^Ci), as is evident from 
Fig. [22H . This is possibly due to the fact that the difference between the full O (a^) results 
and the first order expansion of the resummation, shown in Fig. [22b . is small. 



A.l Weighted recombination in NLO calculations 

NLO Monte Carlo calculations for multi-jet processes are highly CPU intensive. Con- 
sequently, one carries out multiple calculations (runs), spread across many CPUs, and 
averages them so as to get the final result. The correct way of determining the average 
is to weight each run in proportion to its number of events. In practice, however, it is 
common for the distribution of each run to contain one or two bins that are "outliers". 

As mentioned at the beginning of section [2] there is an unknown overaU constant in Yi2{v), which 
relates to the NNLO coefficient of the dijet cross-section. This is irrelevant for us here, since we only use 
the derivative of S2(^^)- 
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Figure 21: Comparisons of logarithms predicted by the resummation with the exact O (a^) 
resuhs (a,b) and the O (a^) results from NLOJET++ for T±^g. Shown for the Tevatron energy 
and cuts, with pn > 200 GeV. 
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Figure 22: Comparisons of logarithms predicted by the resummation with the exact O (a^) 
resuhs (a,b) and the O (a^) results from NLOJET++ for Tm,g. Shown for the Tevatron energy 
and cuts, with pn > 200 GeV. 
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Figure 23: Comparisons of logarithms predicted by the resummation with the exact O (a^) 
results from NLOJET++. The left-hand plot corresponds to the case with (event) number- 
based averaging of NLO results from separate runs. In the right-hand plot, for each bin of 
the NLO results, each run has been given a weight inversely proportional to the square root 
of the error on that bin in the run. The results are for the T±^g observable, with Tevatron 
energy and cuts and pu > 200 GeV. 



obviously inconsistent with the distribution as a whole, and which are a consequence of a 
handful of real and subtracted NLO events with very large opposite-sign weights that end 
up in different bins. These outliers lead to visible anomalies also in the number-weighted 
average and make it almost impossible to use the final distribution without some (often 
questionable) prescription to deal with the outlying bins. 

A common alternative to number-weighted averaging is, for each bin of a run, to choose 
a weight that is inversely proportional to the square of the bin's error in that run. This 
is an option for example in NLOJET++ (and is implicit also for the total cross section in 
programs like MCFM [95] that use VEGAS). Since outlier bins tend to have much larger 
errors than normal bins, they contribute little to the average, resulting in much smoother 
final distributions. However, the error-weighted averaging procedure introduces a bias, 
because there tends to be a correlation between the value in a bin and its error: for 
example, in event samples with positive-definite weights, it is well known that runs with 
larger bin values also have larger errors, and the final error-weighted average systematically 
undershoots the correct result. 

Fig. |23] shows the analogue of fig. |2T1 comparing event number-based and error-based 
weighting. At large negative values of L there is a clear slope, i.e. the bias in the error- 
weighted procedure causes the result to disagree with the expectations based on resumma- 
tion. Only with number-weighted averaging does one obtain results like fig. EH which show 
agreement between the logarithmic structure of the NLO and resummed calculations. 

So as to deal systematically with the issue of outlying bins figs. [2T]and|22]use a modified 
version of the number-based weighting, as follows. One first determines an error-based 
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average for a bin b^'^^ , and a corresponding uncertainty on its contents o"'-"'^ — this provides 
an estimate for the correct value. One then carries out the number-based average with 
the following modification: for a given bin, one excludes runs whose result is further than 
Na^""^ from the fo^"") (we use N = 100 for 15 runs; N should scale as the square-root of the 
number of runs). This then gives us a final result that is smooth and with a substantially 
reduced bias relative to an error-weighted recombination. 

Note that in the phenomenological plots of sections[3]and|l]we have used the error-based 
recombination weights. On one hand the bias that it introduces is modest compared to 
uncertainties from subleading effects. On the other, some of our runs used Rambo [96j phase 
space and others the dipole |97] phase space, and this automatically privileges whichever 
of the two gives best convergence in a given phase-space region. 



B Comment on effect of forward rapidity cut 

For both generic global event shapes and those with an exponentially suppressed forward 
term, in order to satisfy the globalness requirement needed for the NLL resummation, 
we included all particles in the event, including those in the forward/backward regions. 
Experimentally however, it is not possible to perform measurements up to infinite rapidity. 
At the Tevatron the forward detector coverage goes up to ?/ ~ 3.5 and, at the LHC, 
measurements up to y ~ 5 are viable. Theoretical arguments suggest that as long as the 
event-shape's value is not too small, the effect of not including forward emissions should 
be negligible [22], specifically if v > Wmm, with v^^^ given by [23] 

Vmin ~ e-("+^^-2)^-^'' , (B.l) 

where the a and bi parameters were discussed in section 13.1.11 Examining the pure re- 
summed distributions in [23], we came to the conclusion that the result in eq. (IB.1|) for 



fmin ensured that the cutoff would usually have an impact only well below the maximum 
of the distributions. Here we supplement this analysis with a numerical study that inves- 
tigates the impact of the rapidity cut in practice. 

For this purpose we compute the NLO-I-NLL prediction using a rapidity cut on input 
particles for the NLO part of the calculation and compare this to the full NLO+NLL 
without forward rapidity cuts. In parallel we carry out an estimate using a Monte Carlo 
event generator, since it is straightforward to run it with a rapidity cut. 

Fig. EH shows comparisons between NLO+NLL with (solid line) and without the cut 
(full uncertainty band), as well as the corresponding Monte Carlo predictions obtained 
with Herwig (without UE) at parton level at the Tevatron {pn > 200 GeV). Fig. [25] shows 
the corresponding results at LHC {pa > 200 GeV). The results with the rapidity cut 
are always contained in the full uncertainty band of the results without. Furthermore, 
there is in general very little difference between the two Monte Carlo predictions, with the 
exception of the directly global transverse thrust, which is the observable most sensitive 
to forward emissions, as the weight of emissions in the forward region is large compared 
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to that of emissions inside the jets. We note that this is also the one observable where the 
difference between Monte Carlo predictions and NLO+NLL is largest. 
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Figure 24: Resummed distributions for a selection of event shapes with and without forward 
rapidity cut for the high-pj sample at the Tevatron. 
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LHC, 1 4 TeV NLO + NLL (all uncert.) YTTTTm 

Pti > 200 GeV, Ivjetsl < 1 , ric = 1 .5 NLO (|y| < 5) + NLL 

Herwig all partons 

Herwig partons with |y| < 5 

Figure 25: Same as fig. |2l]for the low-pj sample at the LHC. 
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